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Abstract 






We provide dramatic evidence that 'Mellin space' is the natural home for cor- 

Qh relation functions in CFTs with weakly coupled bulk duals. In Mellin space, CFT 

correlators have poles corresponding to an OPE decomposition into 'left' and 'right' 

sub-correlators, in direct analogy with the factorization channels of scattering ampli- 

p^ tudes. In the regime where these correlators can be computed by tree level Witten 

J> diagrams in AdS, we derive an explicit formula for the residues of Mellin amplitudes 

^^ at the corresponding factorization poles, and we use the conformal Casimir to show 

.^ that these amplitudes obey algebraic finite difference equations. By analyzing the re- 

^-H cursive structure of our factorization formula we obtain simple diagrammatic rules for 

f — . the construction of Mellin amplitudes corresponding to tree-level Witten diagrams in 

^^ any bulk scalar theory. We prove the diagrammatic rules using our finite difference 

I equations. Finally, we show that our factorization formula and our diagrammatic rules 

IIV morph into the flat space S-Matrix of the bulk theory, reproducing the usual Feynman 

• ^H rules, when we take the fiat space limit of AdS/CFT. Throughout we emphasize a deep 

rN analogy with the properties of flat space scattering amplitudes in momentum space, 

C^ which suggests that the Mellin amplitude may provide a holographic definition of the 

fiat space S-Matrix. 
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1 Introduction 

Progress in physics often requires a formalism that makes both the symmetries and the 
dynamics manifest and simple. For example, recently we have seen dramatic progress in 
S-Matrix theory facilitated by the use of the spinor-helicity formalism, on-shell superspace, 
and twistor space. We will argue that the Mellin representation [H [21 [3] is the most nat- 
ural framework for CFT correlation functions, especially in the large N expansion. The 
benefits of adopting 'Mellin space' are structurally identical to the crucial, yet rather pedes- 
trian progression from position to momentum space for correlation functions and scattering 
amplitudes in flat spacetime. 

Arguably the most important dynamical property of a conformal fleld theory is its obe- 
dience to the operator product expansion [U [5] [6], [7j , which says that for any two operators 



Oi and O2, we have 

00 

o,{x)om = Y.Y. c^2^Z.^.ix)9,. ■ ■ ■ 9,r.o;--^{o) (1) 

p m=0 

as an exact operator relation, for some coefficient functions Ci2p that are kinematically 
determined up to an overall coefficient for each conformal primary field Op. If we begin with 
an n-operator correlator and use the OPE to reduce it recursively by squeezing together 
k and n — k oi the original operators, then we naturally end up with something akin to a 
factorization channel. 

The Mellin representation displays these factorization channels as poles. We will argue 
that the residues at these poles are intimately related to lower point correlation functions, 
so that the Mellin amplitude inherits a recursive structure from the OPE. We support this 



with a somewhat formal argument in section 2.1 and then develop it extensively for large N 



CFTs in section 3.1 However, the presence of these poles follows simply from the definition 



[H 12] of the conformally invariant Mellin amplitude Mn{6ij) in terms of a CFT correlator via 
{Oi{Xi) ■ ■ ■ On{Xn)) = f[dS]Mn{6^,) J](x, - Xj)-^''^T{S,,) (2) 

-^ i<j 

The integration variables Sij are the CFT analogue of the kinematic invariants Pi ■ pj in 
scattering amplitudes, and we will explain the precise definition of the contour integral in 
the next section. Terms in the OPE of the correlator involving an operator of dimension A 
will have definite power law dependence on the Xi coordinates, and this specific power law 
can only be reproduced by the Mellin amplitude if M„(5jj) has poles on which an appropriate 
linear combination of the 6ij variables can be localized. 

The power of the OPE is that, in principle, with knowledge of the spectrum of operators 
and the three point functions, we could recursively compute all of the correlation functions 
in any CFT. Mellin space may play an interesting role in this program, but we will follow 
an easier route and study correlation functions in large A^ CFTs, specifically those with 
an AdS dual. The Mellin representation becomes far more powerful in CFTs with a large 
A^ expansion and a small number of low-dimension operators [3], because in these theories 
multi-trace operators are equal to products of single trace operators, up to 1/A^ corrections. 
This fact is in a sense built into the Mellin space formalism, as we will explain in section 



2.4, so that the contributions of multi-trace operators are automatically incorporated once 



single-trace exchanges are correctly reproduced. 

The AdS/CFT correspondence [HI El [ID] has led to a long list of major insights into 
both quantum gravity and gauge theory. In particular, it enabled the computation of 4- 



point correlation functions in strongly coupled conformal field theories using supergravity 
[ni[T2l|13l|l2|15l|16lll7l|18l|19l|20l[ai[22l|2g see also ^ for an interesting related 
conjecture. However, progress on holographically computing correlation functions with more 
than four external points has been relatively modest. (Some progress was reported in [25].) 
We hope to facilitate progress in this direction by deriving a precise factorization formula 
that makes it possible to recursively compute tree-level Witten diagrams for arbitrary bulk 
scalar field theories. As a first demonstration of the power of our methods, we compute 
the 5-pt and 6-pt correlation function in 0j0j(/)fe theory in AdS. For example, the Mellin 
representation of the unique 5-pt Witten diagram in 0^ theory with A = d = 4 is 

^'^^''^ ""[jj,^ 36 J6, - 1) + 3(6. - 1)5, ^m.-m,-l)) ^^^ 



where 5a and 5h are linear combinations of 5ij that are analogous to the kinematic invariants 
that appear in familiar fiat-space propagators, such as (pi + P2Y and (p4 +^5)^. 

The fact that our factorization formula can be applied to any combination of factorization 
channels in any order strongly suggests that there must be an even simpler structure underly- 
ing the Mellin amplitudes. In fact, we will derive a set of diagrammatic rules, the AdS/CFT 
equivalent of the momentum space Feynman rules, which permit a direct construction of the 
tree-level Mellin amplitudes in any scalar theory in terms of vertices and propagators. For 
example, in section [5] we define propagators and vertices for a general (pi4>j(f)k theory 

^_^ , VA,A,AA'^i,mj,mk) (4) 

which allow for a direct diagrammatic computation of AdS/CFT correlation functions. Sim- 
ilar results should hold for vector and tensor theories, although aside from a brief discussion 
in Appendix [B] we leave these developments for future work. 

We emphasize that these rules are universal and diagrammatically local, so that the rule 
for one vertex in a given diagram does not depend on the rest of the diagram. At the 
vertices we have 'dimension conservation', which is analogous to momentum conservation 
in fiat space and follows from conformal invariance. The most familiar manifestation of 
'dimension conservation' is the fact that CFT two point functions vanish unless the two 
operators have the same dimension. 

Although OPE factorization may be the physical principle behind the simplicity of 
AdS/CFT computations in Mellin space, this simplicity has another guise in the form of 
a functional equation, which we will derive in section 2.3[ An analogous equation was used 



in position space in |25j and it enabled major progress in the computability of AdS Witten 



diagrams. Furthermore, it was used by [2S] to find explicit expressions for conformal blocks 
with external scalar operators. The main idea is to use the fact that bulk to bulk propa- 
gators in AdS are Green's functions of the conformal Casimir. One can use this Casimir to 
collapse bulk to bulk propagators into AdS delta functions, reducing AdS exchange diagrams 
to contact interactions. In Mellin space, this observation becomes the AdS analogue of the 
very useful fact that V^ — > —p^ when we Fourier transform to momentum space. Namely, 
when we apply the conformal Casimir to the Mellin representation, we find an enormous 
simplification, and a complicated differential equation becomes purely algebraic. 

This functional equation has a variety of uses. We will use it to prove the validity of our 
recursive calculations and as a tool to obtain very general information about the analytic 
structure of Mellin amplitudes. It is also an interesting tool for computing the conformal 
block decomposition of various tensor structures, as we will describe in a forthcoming paper. 

Mellin space has formal advantages because unlike CFT correlators in position space, 
Mellin amplitudes are simple meromorphic functions of their arguments. Again, we have 
a nice analogy with the S-Matrix in flat space, whose analyticity properties are very well- 
known and well-studied. However, on the basis of the existence of a convergent OPE one may 
expect that Mellin amplitudes will always be meromorphic functions without branch cuts 
[1]. In AdS this is reflected by the discrete spectra of quantum theories in AdS, so infinite 
sums replace the phase-space integrals that one encounters in fiat space. The presence of 
poles and the absence of branch cuts was discussed in a 1-loop example in [3] . 

We aim to show that Mellin space is a profoundly useful arena in which to study CFT 
correlation functions, but it may also illuminate the way in which the very large N and 
A limit of AdS /CFT morphs into a holographic description of fiat spacetime [271 [28] , as 
several of us have discussed before |29l |30l [3l [3T|, [321, 133] . We obtain the flat space S-Matrix 
from AdS/CFT by studying high energy bulk states, which are dual to high dimension 
CFT operators, and so we should expect scattering amplitudes to be related to the large 
5ij behavior of the Mellin amplitude. In fact, as was already argued in [3J, we find that 
at large 5ij the Mellin amplitude becomes the scattering amplitude of the AdS theory. In 
section |6] we show explicitly that our AdS/CFT factorization formula reduces to the usual 
factorization of tree- level amplitudes on their propagators when we take the flat space limit. 
This can be viewed as a constructive proof that one can compute the complete tree-level 
fiat space S-Matrix for scalar theories using only CFT correlators, addressing the issues of 
|34| [55| I5B] from a different angle. We also show that the fiat space limit works equally 
well for massless and massive particles. Therefore, we expect that any tree level scattering 
amplitude can be obtained as the flat space limit of the corresponding Mellin amplitude. 



The same should hold at loop level as well, as suggested by the 1-loop example studied in 
[3] . It is then natural to give a holographic and non-perturbative definition of gravitational 
scattering amplitudes as the fiat space limit of the Mellin amplitudes of the dual CFT. 

The outline of this paper is as follows. In section 2, we will use the operator product 
expansion to motivate the Mellin space approach to CFT correlation functions and discuss 
how the Mellin representation becomes especially useful in the large N limit. We will then 
discuss the Mellin representation of a generic scalar Witten diagram and derive the afore- 
mentioned functional equation for the Mellin representation of the diagram. In section 3 
we derive the factorization formula which allows us to split an arbitrary scalar tree-level 
Witten diagram into a 'left' and a 'right' piece. This factorization formula a priori only 
gives us the poles (and residues) of the Mellin amplitude but we will claim that this is in 
fact the entire result. This claim is further motivated in section 4 by working out several 
examples. In section 5 we provide an actual proof of our claim, by demonstrating that it 
is equivalent to the existence of specific Feynman rules in Mellin space and then showing 
that these Feynman rules satisfy the functional equation derived in section 2. In section 6 
we verify that in the fiat-space limit our Mellin amplitudes reproduce scattering amplitudes. 
We end in section 7 with a discussion. Three appendices discuss a more direct proof of the 
factorization formula of section 3, the extension to the exchange of bulk fields with spin and 
some technical developments. 

While this project was being completed we became aware of the interesting work [37], 
which has some overlap with the present paper. 

2 Motivating Mellin Space 

In this section we will explain why Mellin space [IJ makes the physics of CFT correlation 
functions simple and transparent, in the same way that momentum space simplifies scattering 
amplitudes in fiat spacetime. We will see that the key property of scattering amplitudes in 
momentum space is also present for CFT correlation functions in Mellin space: factorization 
of correlation functions on propagator poles, with the residues given in terms of correlation 
functions with fewer operators. Furthermore, in CFTs with a large A^ expansion, the Mellin 
representation of the dual Witten diagrams obeys a simple algebraic equation. We will 



explain these two properties in sections 2.1 and |2.3[ with a brief interlude to introduce 



some notation. Finally in section 2.4 we will review [3J why the Mellin representation is 



particularly well suited to theories with a large N expansion. 



2.1 Factorizaton of CFT Correlation Functions 

First, let us motivate the Mellin transformation by studying one of the most basic properties 
of a CFT, namely the Operator Product Expansion (OPE). 

Consider a CFT correlation function of n operators 

Ik n \ 

An{x,) = lT[0,{x,) J] 0,{x,)\ (5) 

\j=l i=l+k I 

where we have divided the operators into two groups. (We will in this paper only consider 
correlation functions of scalar operators.) Upon recursively applying the OPE, we can write 
the product of operators as a sum 

k oo 

n o, (x,) = Y.Y. ^P,ti-":(^i' • • • ' ^^)^^^ ■ ■ ■ 9,r.o;'-'^{x,) (6) 

i=l p m=0 

where p labels primary operators. 

The general idea behind the OPE is that we can expand in the distance between two 
operators in the limit that this distance is small. Such an expansion is conceivable in any 
quantum field theory, but in a Euclidean CFT we expect the sum to have a finite radius of 
convergence because scale invariance implies that results for small separation should continue 
to hold (as long as there are no other operator insertions). We can apply the OPE to different 
sequential combinations of operators, and the equivalence of these different expansions, i.e. 
crossing symmetry, provides a very powerful general constraint on CFTs [381 ES] • 

Consider now the action of a dilatation on the first k operators, after using translation 
invariance to set x^ = 0, 

k oo 

n O, {e-'x,) =J2J2 e-'(^^+'")+'^' ^'C;-;f:-(a:i, • • • , 0)d,, . . . d,„p;--^{0) (7) 

where Ap is the dimension of the primary operator Op^'"'^^ (0). Since this is an exact operator 
equation (as long as A is large enough to force the first k — 1 points closer to x^ = than 
any other point), we can substitute it into our original CFT correlation function to find 

k 71 \ oo 

l[<^^{e~'x,) n 0.(x.))=5^5^e-^(^^+-)+^S»^»F,,„(xi,...,a;„) (8) 



i«=l i=l+k I p m=0 



where 



F,,^(xi,...,x„) = Cp^-:':-(xi,...,o)(9^,...9^„o;:--Ko) H ^^(^^)) (9) 



i=l+k 



Our question: in what variables is the structure of the OPE manifest? 

A natural answer is to use the variables that are conjugate to the dilatation parameter A. 
In these variables the CFT correlator will have a pole with residue given in terms of lower- 
point correlators. To implement this philosophy, one introduces the Mellin representation 

/n 
[d5]M^{5,,) \[{x, - a;,)-2'5-T(5,,) (10) 

i<j 

where the parameters Sij are symmetric in ij, but (5jj = 0, and they are constrained to give 
the correct behavior under conformal transformations. This means that 

^5., = A. (11) 

i 



Taking into account these constraints, the symbol [dS\ in (10) denotes an integral over a 



subset of precisely n{n — 3)/2 of the 5ij which are independent of each other, normalized as 

The contour of integration for each of the independent 5ij runs parallel to the imaginary 
axis. An extremely useful analogy that will pervade what follows is to think of the 5ij as 
kinematic invariants pi-pj in an n-particle scattering amplitude, and to think of the Aj as the 
masses of these n particles. Then the constraint eq. [TT] follows simply from the requirement 
of momentum conservation ^jPj = and the on-shell conditions pf = — Aj [1]. We will 
discuss below why it is especially natural for theories with a large N expansion to include 
the T{6ij) factor in the definition of the Mellin amplitude. 

Now if we rescale the Xi — ?■ e^^Xi for i < k as above and consider the large A limit of the 
Mellin representation, we find 



i<j i<j i<k<j k<i<j 

To match the leading behaviour at large A between the Mellin amplitude and our OPE result, 
we consider the expansion 

n oo 

n (x|-e-^2x,-x,+e-2^x^)-''^ =^e-^^Q,(xi,...,x„) (13) 

i<k<j q=0 

where Qq is a polynomial of degree q in Xi with i = 1, . . . ,k. Therefore, the contribution of 
a spin / operator to the OPE, comes from the q = I term in this expansion. Matching the 

7 



e ^ scaling with ([s]), we conclude that the Mellin amplitude must have poles at 

k k 

5]]Ai-2 J] 5ij = rp + m (14) 

for all non-negative integers m. Here, we have introduced Tp = Ap — Ip, the twist of the 
operator Op^"'^K Notice that the left hand side is the precise analog of the fiat space kinematic 
invariant —{pi + ...+ PkY- Corresponding poles arise explicitly when we consider Witten 
diagrams in AdS/CFT, and a major goal in what follows will be to give a precise and 
computationally useful formula for the residues of these poles. 

But let us first give an intuitive explanation for why these residues should be intimately 
related to lower point correlation functions. The residue corresponding to a specific OPE 
channel is most conveniently written by introducing for every primary field Op a correspond- 
ing shadow field Oj^ . defined such that: 

{Op{x)dp,{y))=5\x-y)5p,p, (15) 

Clearly, if Op has scaling dimension Ap, the shadow field must have scaling dimension 
d — Ap. An intuitive way to write the shadow field is via the convolution: 

^.(^)=/^'^ (,_'^;/'L.) (16) 

but formally this integral is divergent and needs regularization. 
Using the OPE, we find that at least schematically 

^n(x,) ~5^ /"rf'^y(no.(^i)0p(i/)) loM n ^^(^^)) (17) 

This equation is however only formal as the integral over y of the insertion point of Op{y) 
implies that we may destroy the convergence of the OPE of the other operators. Nevertheless, 
it can be used to offer a reasonable CFT intuition of the OPE in Mellin space. In particular, if 
we were to substitute the Mellin transform of the two correlation functions on the right-hand 



side of (17), the resulting Mellin transform of An has poles precisely at (14). These poles 
isolate specific terms in the sum, and have residues which are given in terms of the product 
of Mellin transforms of the lower point correlators. In section [3] we will see an explicit and 
precise confirmation of this rough OPE intuition in the case of Witten diagrams. 



^Very roughly speaking, one introduces shadow fields in order to write the operator 1 as a sum of 
primary operators acting on the vacuum, ^ Cp|0) (0|Op. Shadows are necessary to ensure that the correlator 
transforms correctly under dilatations; their necessity is analogous to the fact that on a certain very formal 
sense, the bra and ket in-states {p\ and \p) have opposite energy. 



2.2 Conformally Covariant Notation 

We will be discussing CFT correlation functions, so it is natural to use variables [101 SI] 
that are acted on linearly by the Euclidean conformal group 5*0(1, d + 1). If we begin with 
{d + 2)-diniensional Minkowski spacetime, then the conformal generators will simply be 

SO that conformally invariant functions can be constructed out of the covariant inner products 
(Xj ■ Xj). We can view Euclidean AdS(i+i as the hyperboloid 

x^ = -R^ , x° > , X e yf+^ , (19) 

embedded in this {d + 2)-dimensional Minkowski spacetime. We set i? = 1 in what follows. 
Furthermore, we can think of the conformal boundary of AdS as the space of null rays 

p2 = , P ~ AP (A G M) , P e M'^+2 . (20) 

Then, the correlations functions of primary scalar operators of the dual CFT are encoded 
into S0{1, d+1) invariant functions of the external points Pj, transforming homogeneously 
with weights —A,,. We will work extensively with this formalism in what follows; for more 
thorough discussions, see [il] . 

Using the standard AdS/CFT prescription, CFT correlators can be computed in terms of 
Witten diagrams, which are bulk Feynman diagrams that connect to propagators ending on 
the boundary of AdS. The external legs in such diagrams represent AdS bulk to boundary 
propagators, which in this notation are simply given by 

GMP.X) = ^-,^^—^, (21) 

where 

^ r(A) _ T{h + c) 

^ 27rir(A-| + l) 27r^r(c+l) ' ^ ^ 

Here, we have used the notation 

A = /i + c, d = 2h. (23) 

This normalization was obtained by taking the limit of the bulk to bulk propagator; it differs 
by a factor of (2 A — d)~^ from the normalization of |15j . 



Using again the results of [15], this implies that in our conventions the two-point functions 
are normalized as 

(Oa(Pi)Oa(P2)) = ^_J^^.p^^^ ■ (24) 

9 



Finally, to recover the usual expressions in physical M'^ or M'^ we choose the light cone 
section 

P=(P+,p-,P'^) = (l,x^x'^), (25) 

where fi = 0,1, . . . ,d—l. Then P12 = — 2Pi • P2 = (xi —X2Y- We will use Pj and Xi variables 
interchangeably in what follows. 

2.3 The Functional Equation, or AdS/CFT Turns Algebraic 

Let us begin by giving a brief and tortured-looking review of why momentum space drastically 
simplifies tree-level computations of S-Matrix elements in flat space, using 0^ theory as a 
simple example. A directly analogous procedure will lead to a beautiful functional equation 
for the Mellin representation of correlators in AdS/CFT. 

Consider an n-pt correlation function for 0^ theory in fiat space. Isolating the class of 
position-space diagrams with a propagator connecting particles 1 and 2 to the rest, we have 

Thus if we act on this correlator with (Vxi + V^a)^, we can exchange this operator inside 
the integral for V^, which collapses the propagator ^^2(l_ )2 to a delta function. Now if we 
consider the Fourier transform C{pi) of the correlator, this means that 

{pi+P2fCipi) = Co (27) 

where Cq is the same set of Feynman diagrams with the dependence on particles 1 and 2 
reduced to a contact interaction. The transition to momentum space therefore allows us to 
solve for general tree level correlation functions without doing any integrals. 

Now let us perform the analogous steps for correlation functions in AdS, with the Mellin 
representation playing the role of the momentum space amplitude. This procedure has been 
used in position space by [25|. Consider a Witten diagram where particles 1 and 2 are 
connected to the rest of the diagram by a 0^ vertex and a bulk-bulk propagator as shown in 
figure [1| 

AiP,, P2, ..., P.) = ^^^ ^-^^^ (_2P^^X)Ai (-2a\)A2 ^^(^' ^)^^^' ^3' •- Pn) (28) 

where Ga{X, Y) is the bulk to bulk propagator for a field with dimension A. The equivalent 
of the box operator in fiat space is the Casimir of the conformal group, which is just the 
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VLs + A(rf-A) 




Figure 1: By acting with the conformal Casimir on a Witten diagram with a bulk to bulk propaga- 
tor, we collapse the propagator into a delta function. We derive the functional equation by looking 
at this process in Mellin space. 

sum of the squares of the generators. These take an especially simple form in terms of the 



Pi variables, as we saw above in equation (18). The Casimir for the first two particles is 



(29) 



and when it acts on the correlator, inside the integral it is equivalent to ^ J^ = — V^^g. Since 
the bulk to bulk propagator is the Green's function of this operator, 

[VLs - A(A - d)] Ga(X, Y) = -6{X, Y) , (30) 

it collapses Gi^{X,Y) into a delta function. This gives an equation 



■(Ji + J2)'-A(rf-A) 



A = A, 



(31) 



where in Aq the propagator has been collapsed into a contact interaction. In ^25j, this was 
used to convert Witten diagrams with bulk to bulk propagators to contact interactions. 



In Mellin space, this equation takes a remarkably simple form. When the conformal 
Casimir of particles 1 and 2 acts on the product Y\i^{,Pij)^^^^ in the definition (10) of the 
Mellin amplitude, where Pij = — 2Pj • Pj, we find 



{5,n - A)(rf - A - 5,n) + E ^^^'^^^^ (^ " ^ + 9^) 



n(^^^) 



(32) 



i<j 
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where Slr = A1 + A2 — 2(5i2 is the natural analog of the momentum space variable —(^1+^2)^; 
later on we will see that the Mellin amplitude has poles in this Slr. This expression can 
be simplified by noting that multiplication by the kinematic invariants Pij is equivalent to 
shifting the 6ij, so that for example 



*j 



j^^A(p,....,p„) = J m ^ (4,3 - 1)(«. - 1) ) nr(^«)^« 

This allows us to write the functional equation 

{6lr - A)(rf - a - 6lr)M + Y, 2 {Sl^S2jM - d^.d^^MH^ + dy^d^.M^^) = Mq (33) 

where we define 

M^^ = M{Si2 + 1, S^, + 1, Su - 1, S2j - 1, ... ) (34) 

and analogously for the other indexed Ms. This gives a purely algebraic equation for any 
Mellin amplitude M with a propagator connecting particles 1 and 2 to the rest of the diagram. 
Because of the finite differences this equation is more intricate than the analogous equation 
in momentum space, but it will be extremely useful later on for proving general results. In 
particular, in section [5] we will derive a set of Feynman rules for Mellin amplitudes and use 
the functional equation to prove that they correctly compute Mellin space Witten diagrams. 
For completeness, the general functional equation corresponding to any propagator is 

Mo = (Slr - A) (rf - a - Slr)M + Yl {^a^5bjM - 5aj5uM'^^ + (5„fe(5,,M2^) (35) 

ab<.k<ij 

where the propagator separates the first k from the last n — k operators. 

We should also point out that the functional equation is useful beyond its application to 
Witten diagrams. To see this, return to the OPE equation ([6]). Acting with the conformal 
Casimir on both sides, we see that whenever the operator Op on the right belongs to a 
conformal representation of lowest weight Ap and spin Ip, we find: 

^ k k oo 

2 E '^^)' n ^^ (^«) = - E 9i\^ ^v) E ^;;.;--":(^i' • • • ' ^'^)^m • • • d^rr^O^p^-^'M, (36) 

i=l i=l p m=0 

where the conformal Casimir g{A,l) = A(A — d) + l{l + d — 2). So, the contribution of a 
primary Op and all its descendants to a correlation function can be packaged into a single 
solution of the homogeneous functional equation. For the four point function, this solution is 
determined entirely by the kinematics. It was used by Dolan and Osborn to find explicit and 
simple expressions for these contributions, which are the familiar conformal blocks |l2l [26] . 
In Mellin space, conformal blocks look even simpler, as we will describe in a forthcoming 
paper where we will perform a general analysis of the solutions of the functional equation. 

12 



2.4 The Mellin Representation at Large A^ 

We have motivated the Melhn amphtude as a natural representation of CFT correlation 
functions which makes the structure of the OPE manifest and obeys algebraic equations 
whenever there is a local bulk description available. In this section we will explain why we 
have included the factors of T{6ij) in the definition of the amplitude. 

Consider a large A^ CFT with single-trace primary operators Oi and O2 of dimension Ai 
and A2. When we investigate their OPE, we will always have contributions from operators 
such as Oid'^"'02 with dimension Ai + A2 + 2?2 + 0{jj). The product of P functions in our 
definition of the Mellin amplitude then guarantees that the full Mellin integrand has poles 
whenever 6ij = —riij for integers riij. In our simple example, the poles at 612 = —riu produce 
residues of the integral that combine Oi and O2 into an operator of dimension Ai + A2 + 2ni2, 
exactly as desired. 

In general, we can combine k operators in the same way. Taking the first k for conve- 
nience, we find poles at Sij = —Uij for all i,j < k. This set of coincident poles produces the 
multi-trace operator 

where the derivatives are contracted so that between each i and j there are Uij contractions. 
The residue at this pole is related to the correlation function of this operator with the other 
n — k operators via our factorization formula, giving a sort of LSZ prescription for extracting 
the correlation functions of many composite operators. 

Notice that the above argument is very specific to CFTs with a perturbative description. 
In more general CFTs the analogous 'multi-trace' operators have finite anomalous dimensions 
and the poles of the T{5ij) factors are not at the correct location to account for the multi- 
trace contributions to the OPE. In those cases M{5ij) must not only have additional poles 
accounting for the multi-trace operators but also zeroes to cancel off the poles from the 
P((5jj) factors [Ij. In our case the T{5ij) already give poles at the right location and we do 
not need to worry about zeroes in the Mellin amplitude. Perhaps a different definition of 
M{5ij) would be appropriate for such theories. 

Of course multi-trace operators do gain anomalous dimensions beyond the leading order 
in 1/A^, which can be read off from the leading connected higher-point correlation functions. 
The anomalous dimensions of the multi-trace operators can be obtained from the Mellin 
amplitudes. To see how this works, note that a contour integral around a double pole gives 

d5 x-^ 



2Tii {5 - A)2 " 
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—X log X . (38) 



Logarithms are just the perturbative manifestation of the anomalous dimensions, because 
x''' ~ 1 + 7logx. So anomalous dimensions can be extracted by studying the double and 
higher poles of the Mellin amplitude. As a concrete example, the 4-pt amplitude in g{(t)xY 
theory has factors of T{5ij) in the Mellin integrand. So in the OPE channel combining O^ 
and O^ as operators 1 and 2, we have (5i2 = ^34 due to the constraints (recall that this is 
analogous to (pi +P2Y = (ps + ^4)^)- This leads to double poles in the Mellin integrand due 
to r((5i2)r(534), which tell us that the operator O^O^ has received an anomalous dimension 
proportional to g. However, if we look at the double-trace operator O^O^ in a theory where 
A<^ — A^ is non-integral, then there are only single poles, since ^12 = ^34 — A^ + A^. So the 
operator O^ does not receive an anomalous dimension in this theory, as expected. One can 
go on to consider more complicated examples, although it is necessary to disentangle the 
contributions from primaries and descendants. 



The above discussion leads to the general expectation that the inclusion of the T{5i 



Vj 



factors in the defining equation (10) should lead to simple Mellin amplitudes for Witten 
diagrams. This effect is most extreme when we consider the Mellin amplitude of an n-point 
scalar contact interaction in AdS, which is simply a constant |3j. 



^'^ (39) 



I/-4 p^f^ 4<^^-l^) n ill) /i''^! n -(^-^ w. 

This is analogous to the fact that an insertion of a </>" vertex is a constant in the momentum 
space Feynman rules for scattering amplitudes. 

Since this formula will be one of the main ingredients in the remainder of the paper, let 
us quickly review its derivation, which is essentially a slight generalization of a result by 
Symanzik |13]. One first writes the bulk-boundary propagators as 

1 i / U-h Ai-2tiP,-X 



-2P, ■ X)A. r(A,) 7o U 



^p^2un-x (40) 



In this representation, it is easy to integrate over the bulk coordinate X to obtain an expo- 
nential of ^ titjPij. Then the key is to represent some of the terms in the exponential using 
the standard Mellin identity 

|.rw.- (41) 

and evaluate the tj integrals, leading to the Mellin representation with the constraints that 
we have discussed. This is how our standard representation for the Mellin amplitude arises 
in the case of the simplest AdS amplitude. 
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3 Factorization of AdS/CFT Correlators 

In this section we will show how to recursively compute the tree-level Mellin amplitude 
for any theory of scalar fields, including any number of derivative couplings. As modest 
examples, in section |4| we will compute the 5-pt and 6-pt functions in a theory with 3- 
pt contact interactions among scalar fields with arbitrary masses, dual to operators with 
dimensions Aj. 

These computational advances are made possible by a factorization formula that we 
will derive shortly, which says that a Witten diagram with a propagator that divides the 
amplitude into a left and right piece (see Fig. |2]) will have simple poles at 



k 

5lr = a + 2m , where ^^r = ^ A^ - 2 ^ 5ij (42) 

is the direct analogue of the familiar intermediate propagator variable from flat space scat- 
tering amplitudes, — (pi + ■ ■ ■ + PkY- Here A is the dimension and twist of the exchanged 
operator; since we are dealing with bulk scalar fields these are identical. The residues at 

these poles are 

,^T{A-h + l)m\ 



4^ -(A^T^Tltr ^""^^''^ ^""^^''^ ^^^^ 



where 



Lus.,)= Yi <-M(^^.+^«.)ri%V (44) 

Rm{k,)= Yl Mn-k+liS^^ + n^,) U ^-^ (45) 

J2nij=rn k<i<j ■' 

which depend only on the lower point diagrams or Mellin amplitudes, M^^^ and M^]^_^. 
We remind the reader that the Pochhammer symbol (x)„ = T{x + n)/T{x). The derivation 



is given in section 3A The basic idea is to re-write the Mellin representation of a Witten 
diagram with a particular propagator in terms of the amplitudes to the left and right of 
this propagator. Then we massage the propagator into a form such that the entire Witten 
diagram is written in terms of lower-point Mellin amplitudes. This will allow us to identify 
all of the poles and residues in the appropriate 5lr variable. 

Our factorization formula can be applied to any factorization channel. The equivalence 
of all the different possible recursive applications of the formula seems to be a very strong 
constraint on the form of Mellin amplitudes. This suggests that there may exist a single set 
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dc / dP 




Figure 2: A pictorial representation of the derivation of the factorization formula. 



of diagrammatic rules whereby an arbitrary Mellin amplitude can be constructed. In fact, 
in section |5] we will see that such a set of rules exists, and we will prove that it gives the 
same results as the factorization formula. 

Finally, in section [5] we will show that our diagrammatic rules and our factorization 



formula satisfy the functional equation from section 2^, and in appendix |A] we complete the 
proof that our formulas are identical to the Mellin representation of the Witten diagram. 
Thus in any scalar theory, one can compute all tree-level Witten diagrams either recursively, 
or by explicit construction using the diagrammatic rules. 



3.1 Factorization on AdS Propagators 
3.1.1 An Integral Formula 

Our starting point is the the following formula for the bulk-to-bulk scalar propagator [3] 

dc 2c^ /■ ,„ C/j+c Ch-c 



J —I 



dP- 



(46) 



2m c2 - (A - hy J (-2P ■ X)'^+^ (-2P ■ Y)''- 

where we recall that h = d/2. Notice the appearance of two bulk-to-boundary propagators 
in the integrand. Beginning with an n-point Witten diagram, we will use this representation 
for a specific internal propagator. This propagator will break the ra-point Witten diagram 
into two Witten diagrams integrated over a common boundary point, 

"'^ dc 2c2 



An ( J ,; 



27ric2-(A-/i)2 



I dPAl^,{P,,, P)AtuAP^^ P) (47) 



We represented equation (47) pictorially in Fig. pi The left (L) amplitude has k+1 external 
legs with weights A, at point Pj for i = 1, 2, . . . , A; and weight /i + c at point P. The right (R) 
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amplitude has n — k + 1 external legs with weights Aj at point Pi for i = k + l,k + 2, . . . ,n 
and weight /i — c at point P. In other words, the operator with coordinate P on A^ is the 
shadow of the operator with coordinate P on A^. Now we write the left and right amplitudes 
in the Mellin representation 

A^+i(p„p) = Arf5]i[rf/]iM,^+i(4.,/,)nr(4)(-2P,-p,)-^'^nr(/,)(-2P.-P)-'» (48) 
■^ i<j i=l 

/n n 

[MUdl]RMt,+^{l,,U) n r(4)(-2P,-p,)-^-^ n r(/.)(-2P,-p)-'« 
k<i<j i=k+l 

The integration measure [(i(5]L[(i/]i is constrained by 



k 

y^ h = h + c 
1=1 

k 
h + J2^iJ = ^i (^ = 1,2,...,/C) 



(49) 



and analogously for [(i(5]/j[(iZ]/j we have the constraints 

n 

y li = h — c 



i=k+l 



(50) 



/j + ^ % = Aj (i = A; + 1, . . . , n) . 

k<jj^i 



The equations in (48 ) may appear complex, but they encode very simple information, namely 
that A^ is a CFT correlator between the k operators on the left of the propagator and a 
new operator with dimension /i + c, and equivalently for A^. The constraint equations are 
identical to what we would find from momentum conservation if we introduced fictitious pi 
with 5ij = Pi ■ pj and p] = — Aj. 

Inserting now the Mellin representations (48) in (47) we find 

_ 2^^ c^-iA-hy J [^'^]^[^^]^ ^'^'^^'^' ^'^ n r(5^.)(-2P. ■ P,)-'- 

j [d^]R[dl]RMtk+l{^^vli) n r(4)(-2P* ■ P^'-' J c^Pnr(/^)(-2P, ■ P)-'^ (51) 



k<i<j i=l 



Notice that the contour integral over c requires knowledge of the Mellin amplitudes M^_^^ 
and M^_i.^^ for general (complex) external scaling dimensions Aj. As we will see explicitly 
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in all the examples below, for general Witten diagrams we obtain a Mellin amplitude which 
is an analytic function of the Aj and so this is not problematic. 

The last integral was studied by Symanzik |13] and has a simple Mellin representation 

/n „ n 

d'PllT{k)i-2P, ■ Pr^ = n^ / [d6] llTi6,,)Pr/^^ (52) 

where the measure [d6] is constrained by J2]=/=i ^ij = ^«- We have reduced the dependence 
on all external kinematic invariants to the Mellin form! This means that we can shift the 
integration variables 5ij — )■ 6ij — 6ij to obtain the following expression for the Mellin amplitude 
of the original n-point diagram 



Mn{5.^) = / 7^7^^^^^^:^ LxR (53) 



27ii c2 - (A - hy 



where 



k<i<j '■' 

It is crucial that the factor L depends only on 5ij with 1 < i,j < k and the R factor depends 



only on 6ij with k < i,j < n. Using the constraints (50) we can solve for and eliminate the 
li variables, which are the last vestige of the spacetime version of the internal propagator. 
This leaves the single constraint 

k ^ k 

2j2^ij + h±c-J2^i = ^ (55) 

i<j i=l 

which expresses the dimension h ± c oi the internally propagating operator in terms of the 
other integration variables in L and R, respectively. 

Let us end the presentation of the factorization formula with some comments regarding 



the relation with the OPE decomposition as presented in equation (17). At first sight, 
the above factorization is very similar to the factorization appearing in equation (17). The 
difference lies in the additional integral over c in the evaluation of the Witten diagram, which 
in the analogous CFT computation is localized at c = A — h. Although we will not prove this 
statement here, we claim that this additional integral is eventually responsible for correctly 
taking into account the contributions of multi-trace operators. (At the level of the four- 
point function this follows from equation (38) of [3], where the integral over c is explicitly 



responsible for the multi-trace poles in 612 and 634.) Furthermore, because the choice of 
contour for the c-integral breaks the symmetry between the field and the shadow field our 
expressions can be made completely exact and are not of the formal form as in equation 



(17). It would be interesting to make the relation between the multi-trace contributions 
and the c integral precise, for which an analogous discussion of the conformal partial wave 
decomposition in [1] should also be very useful. 



3.1.2 Identifying the Poles 



The advantage of equation (53) is that we can use it to determine the residues of the poles 



of the n-point Mellin amplitude in the factorization variable 

k n ^ f ^ \ 

i=l j=k+l i=l \ jf^i / 

Notice that the equivalence between the two expressions on the right-hand side follows from 



the constraints (11). As we discussed above, 6lr is the direct analogue of the kinematic 
invariant that vanishes at the factorization channel in a scattering amplitude, such as — (pi + 
■ ■ ■ + PkY in an n-pt scattering amplitude. 

Contour integrals have singularities when a pair of poles in the integrand collides and 
squeezes the contour of integration between them. We have a many-dimensional contour 
integral, so identifying such occurrences would naively be a daunting proposition, but our 
task is greatly simplified by the structure of the integrand, its contour, and the fact that 
we are focusing on Slr- In particular, since 6lr involves all of the left or all of the right 
6ij, we can obtain a pole in this variable only if we use the poles from all the T{6ij — 6ij) in 
(54), which form a semi- infinite sequence towards positive real infinity in the 6ij integration 



variables. We obtain a pole in S^fi when these singularities collide with the 'propagator 
singularity' at c = ±(A — h), where c is a function of the 6ij from the constraints. 



Before studying the general case, let us analyze the case where M and M are inde- 



pendent of Sij so we can perform the integrals in (54) explicitly. This case corresponds to 



a Witten diagram with a single bulk-to-bulk propagator connecting a (/c + l)-vertex to a 



(n — A; + l)-vertex. To compute the integrals in (54), one starts by eliminating S12 using the 
constraint (55). Then, all other 6ij are independent integration variables and, for constant 



M^, the integrals can be computed iteratively using Barnes' lemma: 

— / T{a + s)T{b + s)T{c - s)T{d - s)ds = -^ ' ) '' ( ^ . (57) 

-^TT? j_joo L [a + + c + d) 
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This gives 



r(^H^^^)r( 



H+c-Slr 



) 



-M' 



which has poles at c = 5lr — h — 2m, for m = 0, 1, 2, ... , with residue 



(58) 



ml 



r(E 



i<j ^ij/ 



-M' 



(59) 



In this expression, M^ depends on the integration variable c. It is convenient to make this 



dependence explicit. Using (|39|) we find that 



Y (" EtiA^-fe+c 



M'^(c) 



n 



r(A-/i + i)r 



2 



4r(c+i) f^^r(A,) r(i + c)r(^-^Y^'^' 



<+i 



(60) 



where M^^^^ stands for the left Mellin amplitude with physical external dimensions Aj and 
A. Similarly, the right part of the diagram gives 



R 



r (^ ^">fe^'-^+" ^ r (^^^^^) T{A-h + i)r (^ ^">fe^'-^- 



P / Er>fc A'~'^j^fl \ 



r(i-c)r 



Er>fe A,+A-2/. 



•M, 



n-fc+1 ' 



(61) 



which has poles at — c = (5/,^ — h — 2m for ?7i = 0, 1, 2, . . . . The poles in Slr of the n-point 



Mellin amplitude arise from pinching the c integration contour in (53) between two poles of 



the integrand. From Fig. [3} it is clear that, for Slr = A + 2m, the contour will be pinched 
in two places, aX c = A — h and at c = h — A. The contribution to the residue of M„ at 
6lr = a + 2m from the pinching at c = A — h involves the product of the residue of L times 
R evaluated at the pole. This gives a residue of the form 



Ti^{A-h) 2- 



-ly 



m\ 



r(Li<,'^ii) 



5ifl=A+2m 



2n^V{A-h + l)^^, ..ij 



m!(A-/i + l) 



-K+iM„«- 



fc+i 



E^4JE^^ 



i<j 



k<i<j 



(62) 



(5LH=A+2m 



which is symmetric between left and right. The other contribution is equal and just doubles 



this residue. In (62), we are always evaluating the residue at the pole Slr = A + 2m, which 
using (56) implies a constraint for the Sij appearing in (62). Notice however that those 6ij, 



being the arguments of the full Mellin amplitude (53), are in principle not subject to this 



extra constraint and we therefore added the square brackets and the explicit 6lr = A + 2m 
to indicate this extra constraint. 



20 



h-A 


h - (5lr + 2m 
A-h 


-h + Jlr -2m 


li 



Figure 3: Structure of poles of the integrand of the factorization equation (53) in the c 
complex plane. When 6lr = A + 2m, the integration contour is pinched between poles at 
two places. 



Let us now consider the general case. First consider the contribution from the collision 
of poles of L with the pole at c = A — /i in (53). The relevant poles from the left 6ij arise 

from the pinching of 

2tt''c^ 
6ii) and — — rnr (63) 



where we remind the reader that c is linked to the 5ij by constraints. This happens when 
the 6ij variables are such that 



This gives the following contribution to the residue of M„ at the pole 5lr = A + 2m, 



(64) 



2ti^{A - h) 



2_j Tlij — ' 



M,^+i(,5,,+n,,)n 



;-l)"-r(,5i,- + n, 



i<j 



n. 



ly 



r(<5i 



X R 



(65) 



J <5ij,=A+2m 

where it is important that the sum of the Uij are constrained to equal m. Notice that 
evaluating the 5ij at the pole 5lr = A + 2m also ensures that the arguments of Mj^j^^ always 
satisfy the appropriate constraints (which follows from (56)) where M^^i is well-defined. 



'fc+i 



This is almost our factorization formula eq. (43), except that here the arguments in R are 



constrained by eq. (55) with c = A — h. Thus, it has the appearance of a lower-point Mellin 



amplitude where one of the legs has been replaced with a shadow field, exactly as we should 



expect from our discussion in section 2.1 In Appendix CI we prove an identity relating 
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the Mellin amplitude M with a shadow field replacement to the original Mellin amplitude: 






k<i<j 



nk 



5LR=A+2m 



r(A-/i)(-l)™m! 
(A-/1 + I), 



y^ M„_fc+i((5jj +nij) Y\ 



— r ['Jijjnij 



Y^nij=m 



k<i<j •' 



(66) 



J <5£_R=A+2m 



Notice that the arguments of M^_f^,^ again satisfy the required constraints. Inserting this 



identity into equation (65), we obtain half of our factorization formula (43). The other half 



comes from the collision of poles in R with the pole at c = /i — A in (53). 



3.1.3 The Complete Factorization Formula and Its Interpretation 

We have shown that any Witten diagram will have a Mellin representation with the above 
poles and residues in the 5lr channel. If the Mellin amplitude vanishes for large 5lr then it 
would be completely determined by its poles and residues, and we would be able to write: 



M 



E 

m=0 



Res(m) 



>LR 



A -2m 



with 



Res(m) 



47r'^r2(A-/i+l)m! 
r(A-/i + l + m) 



[Lm{5ij)Rm{Sij)] 



SLR=A+2m 



(67) 



(68) 



where L^ and Rm are given in (44). In fact, we will see in all examples that a stronger 



statement is true. Our formula is equivalent to its projection onto all of its poles, not just the 
specific 6lr singularity in the factorization formula, so that all of the explicit Pochhammer 
symbols {Sij)ni can be evaluated at poles. If M vanishes as any propagator goes to infinity, 
then this follows from the simple fact from complex analysis that 



E 



/.(') 



z — a,- 



^—^ z — as 



(69) 



when the sum vanishes as 2; — )■ 00. In what follows, when we refer to our factorization 



formula we will almost always be referring to equation (68) with all 6ij in the numerator 



projected onto poles, because it is this pole-projected formula that we will be able to prove. 
In the remainder of this paper we provide strong evidence that the Mellin amplitude is 



in fact completely determined by its poles and therefore (67) is the full answer. This we will 



do as follows. We will first show that our factorization formula implies a set of diagrammatic 
rules for the computation of Mellin amplitudes, and then we will show that these rules satisfy 
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the functional equation from section 2.3 Assuming that Witten diagrams are polynomially 



bounded at large 6ij, as we discuss in appendix Kl this leads to a proof of the factorization 



formula (67) as well as our diagrammatic rules. Additionally, we provide a more direct proof 



of (67) and a few other technical details in appendix [Aj 

3.2 Adding Derivative Interactions 

The results above generalize to scalar theories with arbitrary derivative interactions, due to 
a beautiful interplay between the interaction vertices and our factorization formula. In fact, 
we will see that all dependence on the 6ij from derivatives passes through our factorization 
formula and simply leads to an overall factor multiplying the Mellin amplitude and the shift 
of a few constant factors. 

As shown in [3], the Mellin amplitude for a contact interaction with an arbitrary number 
of derivatives is a polynomial. A convenient basis for these polynomials is 

k 



M, = gl[{6,,)a,, (70) 



i<j 



where the aij are integers related to the number of derivatives coupling field i to field j. 
(Notice however that the derivation below goes through for non-integral aij as well.) If we 
plug this amplitude into our factorization formula, we find 



Y,nij=m i<j ' *•'' ij 



and equivalently for Lm- In other words, M^ has passed through the Pochhammer symbols 
from the factorization formula, and so we can evaluate it at the original 5ij and not at 
^ij + nij. The only effect of the derivative interaction is to shift the Pochhammer symbols 
in the factorization formula, but since these 5ij are naturally evaluated at poles, this only 
shifts certain constants in the residues at those poles. 

This effect persists under recursion, so Mellin amplitudes with many derivative interac- 
tions are simply given by the equivalent amplitude without derivatives with some constant 
shifts and an overall factor from all of the various polynomial Mellin amplitudes from the 
derivative interactions. In particular, this means that if Mellin amplitudes without derivative 
interactions are entirely determined by their poles, then our factorization formula applies 
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Figure 4: Four-point and five-point Witten diagrams in cubic scalar theory. 

to all scalar theories in AdS. Another way of saying this is that when we add derivative 
interactions, the 'skeleton diagrams' with only the propagators are basically just 'dressed' 
by a polynomial coming from the derivatives at vertices. 

4 Sample Computations 

In this section we will demonstrate the power of our formalism by computing the 5-pt and 
6-pt amplitudes in a scalar field theory with 3-pt interaction vertices. Notice that, as will 
become clear below, using the factorization formula it is even easier to compute amplitudes 
in theories with general V"0'' vertices, since the greatest complication arises from having 
many bulk to bulk propagators. 

Before moving on to a non-trivial computation, let us see how our formalism works in the 
simplest case, that of the 4-pt function. Suppose specifically that we have the bulk interaction 
vertices X4>i4>2(p5 and A030405, and we want the Mellin amplitude for {O1O2O3O4) from 05 



exchange, as shown in Fig. |4[ Applying equation (|67|), we find the Mellin amplitude is 



Mi{5ij) 



E 

m 

E 



-47r''r(A5 -h + l)m! 



5lr - As - 2m {^^ - h + 1) 
1 -27r'T(A5 - /i + l)m 



125" 



m\ 



X 



345" 



{S: 



6, 



m 



(As-Zi + l) 



A 



125 



m!(A 



12,5 j 



34 jm 

ml 

1 



345 



m!(A 



34,5j 



5LR=A+2m 

(73) 



where Aj^ ^ = — — ^^ — - and Aj^-fc is the 3-pt Mellin amplitude for a contact Witten diagram 
with external dimensions A,j, Aj and A^. In the second line, we have used the fact that 
26i2 = -6lr + Ai + A2, 2(534 = -6lr + A3 + A4, and the identity (a - m)^ = jafz;;;- We 
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have also introduced the notation 

2^5 = 6lr - A5. (74) 

This will be convenient in amplitudes with multiple propagators because in those cases there 
will be poles in many different specific linear combinations of the 6ij. In this respect, the 
(5j's are analogous to the combinations (pi + ■ ■ ■ + Pk)'^ + M^ in a flat-space diagram; in 
fact, in terms of the equivalent fictitious Mellin "momenta" pi, 26^ = —{pi + ^2)^ — ^5 = 
~{P3 + ^4)^ ~ ^5- The major difference is that, while fiat-space amplitudes have a single 



pole for each propagator, Mellin amplitudes have the full tower of poles in eq. (67) 



S-Matrix: ^ i — > Mellm: nr^ ^ — , m = 0, 1, . . . (75) 

Note also that the vertices, while more complicated than the simple fiat-space factor 
A, do not actually introduce additional (5jj-dependence, despite the naive appearance of 



equation (68), because (5i2 and ^34 are completely fixed by the constraint 6lb. = A5 + 2m. 
So, continuing the comparison of the four-point function with that of flat space, the vertices 
differ by 

S-Matrix: A^ — > Mellin: '^'^ , m = 0,l,... (76) 

ITT'- \^ij,k) —m 

We will next turn to the evaluation of some higher-point amplitudes, where we will see that 
much of this structure continues to hold. Thus although in section [3] we were only able to 
prove that our factorization formula reproduced the correct poles and residues, in this and 
all other examples that we have computed, miraculous identities and cancellations seem to 
guarantee that this is the entire result. In other words, even the naive version of our formula 



(67), with in particular the explicit Pochhammer factors (5ij)„^ in equation (68), actually 



does vanish at infinity. 

The components of the 4-point calculation will be so ubiquitous that it is useful to 
introduce notation for the factors that accompany the propagator, and for the vertices: 

2»r'T(A-ft + l)m! A.,fc ,„, 

^^<'"> = - (A-ft+l).„ • ''-■-'-'<"' = ™.(A„,)_„.- '^^) 

This is our first hint of the Mellin space diagrammatic, which we will develop in section [5j 
Note that the vertex is not symmetric, because only the leg k is 'off-shell'. We will obtain a 



general expression for this vertex with all legs 'off-shell' in equation (88). 
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4.1 5-pt Amplitude 

The next simplest amplitude is the five-point diagram shown in Fig. |4j We can choose to 



apply our factorization formula, equation (67), on either propagator. In this example, the 



diagram is symmetric so the computation is identical either way, so let us decompose on the 
internal line 7. Then, the left Mellin amplitude is the four-point diagram from the previous 
subsection, which depends only on the Sij through the combination 



26, 



A. 



-26 



12 



Ai + A2 - Ae 



{n 



which is essentially the "momentum" fiowing through line 6. Applying (67) with 267 
6lr — A7 thus gives 

1 -27r^r(A7 - /i + l)m! 



M,i6e,6r) = Y, 



6f — m (A7 



X 



Y^ M4((56-ni2) 



h + l)m 

{612)^2 ('^13)ni3 ('^23)n23 



^12! 



nis! 



^23! 




457" 



{5. 






(.79) 



&'!=m 



The fact that this entire amplitude can be written as just a function of 6, and 6-j is not yet 
manifest, since naively the Pochhammer symbols contribute dependence on 5i2,5i3 and ^23 
separately. However, M4 depends only on 6, = —612 + Ai2,6, so we can first do the sum over 
77-13 and n23 with their sum n = rii^ + ^23 fixed. This is aided by the general identitjjj 



n\ 



E n 



n\ 



0) 



Therefore, only the combinations 612 and 61^ + ^23 actually appear in the 5-point Mellin 
amplitude. But, it is easy to see that these combinations are completely fixed by 6, and ^7! 
The most immediate way to see this is by considering the "momenta" fiowing through the 
diagram: (pi +^2) fiows through propagator 6, and (pi +P2 +P3) = — (P4 +P5) fiows through 
propagator 7, so we have that 



-^45 



-^7 + A45J, 6 



12 



-56 + A 



12,6, 



612 + (^13 + 6; 



23 



-(^e + A 



123,7- 



Now, let us substitute our expansion for the four-point Mellin amplitude in eq. (79): 
1 -27r'^r(A7 - /i + l)m! 



M^{6,,6,) = Y, 



A 



457 



^1^ 



(821 



n\2+n=m k 



^67-171 (Ar-h + l)^ V"^K^45,7)-»n. 

^AiAaAsl^ - ^12)5'a6(^ " '^12)Va3A7A6(^ - ^12) (Ai2,6 - ^e)™^ (A36,7 + 6,- m)n 



6r- k 



ni2'- 



n\ 



^ This identity is easily proven by first noting that X]^o 1 = (1 ^ ^) ^ i ^Lnd then matching powers 



of z in the product ni<i(l + 



{l + z)-^'<^' 
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Figure 5: Left (Right): Six-point linear (star) Witten diagram in cubic scalar theory. 

We have checked numerically that this is a valid, explicit and symmetric formula for the five- 
point function, but we can simplify it by 'projecting it onto its poles', i.e. replace 6q ^ k in 
the numerator. In fact, in section |5] we will explicitly prove that our factorization formula 
holds after this simplifying projection is performed, although, somewhat miraculously, in all 
examples this step has not actually been necessary. 

The projection eliminates the 6ij dependence from the Pochhammer symbols in the nu- 
merator, and we find 



M,{6,,6,) = ^ rA.-Ma._^^ ' 



m,k 



{6q - k){67 - m) 
where we have defined a generalization of the vertex function to include two indicesjj 

^A3A6A7(^,"^) = 



X 



367 



-3-f^2 



-k, —m, -^ 



A3+A6+A7 



h 



(83) 



4) 



A;!(A37,6)-fc"^K^36,7)-m''^ " V ^367 " m, Agyg - k 

This vertex is manifestly symmetric in (Ag, k) 4-7- (A7,r7i), as is required by reflection sym- 
metry of the diagram. Furthermore, one may easily see that it reduces to our earlier vertex 
function Va^a Afc("^) when one of the indices is set to zero: 

^A.A,A,(0,m) = \/A.A,A,(m). (85) 



4.2 6-pt Amplitude 

For the next example, consider the "star" 6-point diagram in Fig. [5] We will apply the 
factorization formula to line 9, so the left amplitude is now the 5-point Mellin amplitude 



^Arriving at this expression for VAsAeArik, m) requires the use of a hypergeometric transformation iden- 



tity. 
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M5{S-/,6s) from the previous subsection. We will see that the final six-point amplitude 
depends only on the combinations 6j, ^g, and ^g, which satisfy 

Si2 = -Sj + Ai2,7, (^34 = -Ss + A34,8, ^56 = Sg + ^55,9 

(^12 + (^13 + ^14 + <^23 + ^24 + (^34 = -^9 + Ai234,9- (86) 

These relations are most intuitively understood by noting that the "momentum" fiowing 
through lines 7,8, and 9 are (pi + ^2), (ps + Pa)-, and -(pi + P2 + Ps + Pi) = (Pb + Pe), 
respectively. 



In evaluating (43), we can again take advantage of the fact that M5 depends only on 6y 
and (^g, rather than all possible Si/s, to immediately sum over ni3,n23,Tii4,n24, subject to 
the constraint ^13 + 7223 + rin + 77,24 = n. The 6-point Mellin amplitude can then be written 

SAs{m) f A56 



ET,j /. . ^ (A12.7 - '^7)ni2 (A34,8 - (^8)7134 (A78,9 + ^7 + (^8 - "^)h 

ni2l riul n\ 

\nr2+n34+n=in 



Upon substituting the expression (83) for M5, we obtain an explicit expression for Mg. 
However, as before, we can simplify further by evaluating 67 and ^8 at the residues of all 
poles This step simplifies the calculation, and it means that we are explicitly using the version 
of our factorization formula that we wrote down in section [3j This diagram is particularly 
interesting because it is the lowest-point amplitude that contains a cubic vertex that connects 
three internal or 'off-shell' propagators. As such, it is the simplest diagram that one should 
calculate in order to obtain the generalization of Va^a a^ to three different, non-zero indices. 
A short computation shows that the Mellin amplitude Mq for this 'star' Witten diagram is 

{l)VAsA4As{k)SAs{k)VA,, 

{6,-l){6s-k){6,-m) 



nj f, X X \ 'S^ VA^A2A7il)SA^il)VAsA4Asik)SAsik)VA,AeA,im)SAgim) 

Me{67,6s,6g) = }_^ _ _ _ VArAsAAl,k,m) 



Lk,m 



where 



TA n 1 \ ^789 

VArA8AA^,k,m) 



X 



E 



/!A;!m!(A78,9 + / + A;)_^(A89,7)-z(A79,8)-fe 

( 'jni2+n( ^)n34+n[ ''^Jni2+n34 



^ ni2!n34!n! 

(1 ~ A78,9 — / — A;)„^2+„34(A89^7 — /)„^2+„(A793 — k)n:^^^n 

Although it is not obvious from the way this expression is written, it is easy to verify 
numerically that it is symmetric under interchange of any two (A, i) dimension- index pairs, 
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e.g. (Ay,/) ^ {Ag,m). Moreover, when any of the indices vanishes, it reduces to the vertex 
functions we have already encountered in the 5-point amphtude: 

^A,A,A,(0,A;,m) = VA,A,A,{k,m). (89) 

So far, the above formula for V^iA Ak ih k, to) is simply another way of packaging the 6-point 
amplitude, and one could reasonable expect that a new vertex function would appear at 
each vertex as one considered higher and higher n-point amplitudes. Surprisingly, this turns 
out not to be the case: the functions SAijn) and Va^a AS^)k,m) are all that is needed in 
order to write down the most general n-point amplitude in 0^ theory. These remarkable 
diagrammatic rules for Mellin amplitudes are the subject of section |5j to which we now 
turn. 



5 Mellin Space Diagrammatic Rules 



In section |3.1| we derived a factorization formula for Mellin amplitudes that can be applied 
to any factorization channel, and in section |4] we used this formula to recursively compute 
several examples. The fact that our factorization formula gives equivalent results when the 
recursive steps are applied in different orders suggests that there exist universal diagrammatic 
rules that allow for the construction of any Mellin amplitude. The purpose of this section is 
to derive and prove these rules. We also obtain a practical benefit, because the computation 
of complicated diagrams becomes standard and straightforward. 

5.1 The Diagrammatic Rules and Factorization 

To begin, we will demonstrate that there is a special case where Mellin amplitudes may be 
computed with essentially exactly the usual fiat space procedure. This case is the A0^ theory, 
in any number of dimensions up to 6y when the CFT operator dual to has dimension A = 2. 
The reason that the usual kind of diagrammatic rules automatically apply to this case is 



that every time we use (43) to add on a cubic vertex to a lower-point diagram, we encounter 
the factor (in, say, R): 

m!(l)_„ 
where M3 is the 3-pt Mellin amplitude. Thus, the sum over m always truncates at m = 0, 
and therefore all the sums over the n^j's do as well. Since all amplitudes in this theory can 



* Above 6 dimensions, A = 2 is not allowed in a unitary CFT. 
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be built up by adding on three-point vertices, there are never any sums to do at all! For an 
arbitrary diagram, the rules for the Mellin amplitude are simply to include factors of: 

1. ^M for each propagator, (where 5i is the appropriate linear combination of 5i/s), 

2. M3 for each vertex. 

So, in this special case, the calculation of Mellin amplitudes is identical to the calculation 
of the corresponding S-matrix elements in flat space, as long as we replace the kinematic 
invariants Pi ■ pj — )■ 6ij. For instance, we can write the five-point diagram from figure |4] in 
terms of the 'Mellin momenta' as 

M,(p„...,p,) = M3 . ~'%'l. M. -'%'l. M,. (91) 

For general A, however, the sum on m will not truncate at m = 0, and each propagator 
will be associated with multiple poles: 

M = T Mim, m.) 

^ ((5i-m,)...(5, -m.) 

where 

S. - -^ (93) 

is the square of the total 'momentum' flowing through the i-th propagator. Thus, the best 
that we could hope for generally is to have diagrammatic rules for the residues M (mi, . . . , m^). 



We have already seen in our examples in the previous section in equations (83) and (88) what 
form such rules might take. Thus we will optimistically guess that more generally, one can 
calculate M (mi, . . . , m^) according to the following rules: 

1- %r^ for each propagator, where SA^^a) = "^''(1^-^+1^"'°' 
2. V'A,A,Afc("T'i,"T'j,'"^fc) for each vertex. 

These rules are depicted schematically in Fig. (pfjj We will show below that this very 
simple procedure is exactly correct for higher-point diagrams as well! On the one hand, this 



is rather surprising from the point of view of the factorization formula (43), where naively 
any diagrammatic vertex factors like VAaA^Ac would have to depend on all the indices of 
all internal lines in the diagram. On the other hand, the existence of such a set of rules is 
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Figure 6: Diagrammatic Rules. The expression for Vabc is given in eq. (88), and that for Sa 



is given in eq. (77). 



very natural in that it automatically explains why the factorization formula gives the same 
answer when applied to any propagator in diagram, a very strong consistency condition. 

The first step is to simplify the Pochhammer symbols. Since all dependence on the SijS 



in Ml and Mr is through the propagator variables, we can always use eq. (80) like we 



did in the case of the five-point function in section [4j There we grouped all the SijS into 
terms that depend on only the propagator variables, the (5jS. However, each 6ij appears in a 
Pochhammer symbol exactly once, so each 6ij can only appear in a single linear combination 
of the 6i. For instance, in the five-point function, we had to use the latter two identities in 



eq. (80) to write S12 and ^23 + (^13 as 



6 



12 



-Se + A 



12,6) 



(^IS + ^23 = -5? + ^6 + A 



36,7- 



(94) 



In general, this regrouping can always be performed, so that for every vertex we have a 
Pochhammer symbol of the form 

{A.ab,c + Sa + Sb — Sc)n^-na-nt 



[Ur. 



Tin 



nb)\ 



(95) 



as depicted in Fig. [7j Here, a and b are the two propagators leading into the vertex and 
c is the propagator leading out, as we work from the external lines inward toward the 
6lr propagator. For completeness, we note that this identity generalizes to arbitrary n-pt 
interactions, with the a and b indices replaced by a sum over all the propagators flowing 
into the vertex, towards the propagator on which we are factorizing. Each Sij shows up in 



^All external lines are taken to have index m = 0. To save space, we will abbreviate S'a„('ti) to Sa{m), 
and similarly for V. 
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(ric-na^ rian)^- 

Figure 7: When evaluating the factorization formula on its poles, products of Pochhammer symbols 
reduce generally to a simple expression associated with each vertex. 

exactly one such Pochhammer, because starting with external vertices i and j and flowing 
inward through the diagram, there is always a unique vertex where the two meet. The sum 
is then reduced from a sum over all the n^j's to a sum over rij's, one for each propagator, 
defined by the condition that 6ij -^ 6ij + n^ is equivalent to 6i — )■ 5i — njjj An important 
nicety of defining the n^'s this way is that the poles are always shifted by 

^ ^^ (96) 



Si — nii 6i — rrii — rii 

As we will soon see, this makes it possible for us to absorb all of the njS via a re-definition 
of the sums. 

To demonstrate how the Feynman rules follow from the factorization formula, let us 
factorize a Witten diagram along a certain 'central' propagator and assume that the asso- 
ciated Ml and M^ have been computed from the diagrammatic rules. We will prove that 



if we compute use the factorization formula (67) to compute M„ from Ml and Mr, then 



the result will be equal to what we would have found had we computed it entirely from the 
diagrammatic rules. 

It is helpful to illustrate our arguments first with an example, so let us consider the 
eight-point diagram below, as it will exhibit all the features we will use for the general proof. 
Factorizing on the "e" propagator. Ml is a seven-point diagram and M^ is a three-point 
diagram. Assuming that we have applied the diagrammatic rules to Ml and Mr, we obtain 
for Ms 



^To be a little more explicit, if pi = (fci + k2 + ■ ■ ■ + ks), where ki are all external "momenta", then 

s 



"i = Z]«^. "' 
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^8 = E 



5e - rUe 



E 



/ 






E 



^12a ("^a) <S'a (ma) V'34b (m^) S'b (mb) V;fec("^a , "^fe , ^c) Vsed (m^) S'd (m^) Krfe ("^e 



{5a -rria- na){6b - rrih - nh){6c -m^- nc){6d - m^ - na) 

(Al2,a — ^a)na (^34,6 ^ ^b)^^ {^ab,c + ^a + ^b — ^cjuc-na-nt {^56,d — ^d)na 



{Acd,e + 6c + 6d- me)„ 



{nc-na-nb)\ 



rid^. 



Ulp 



(97) 



We can eliminate Ug through ric + n^ + rie = mg, and then the sum on nj's is unrestricted. 
To simphfy further, we first redefine rrii —^ rrii — rii in the sums on m^ in order to shift the 
poles back to rrij, and then as usual we evaluate all the 5iS in the numerator on the poles 
(i.e. 5i -)■ rrii). 

Now we want to show that Mg also satisfies the diagrammatic rules. We will first show 
that the correct V and S factors are associated with the residues of the poles in 5a and 5^. 
This follows from the following identity, which we have verified numerically: 

E( M f \c I N (Al2,a - ma)„^\ ( ^. j .^j (A^ib- 'mb] 

I Vi2a{ma - na)Sa{ma - Ua) ^^-j ) ( V^AA^h " nb)Sb{mb - Ub) 



Ual 



rib 



nb\ 



X Vab,c{rna -na,mb- Ub, m^ - nj 



{Aab,c + ma + mb- mc)nc-na- 



"6 



Iric - Ua - Ub l 



(9^ 



T. / xcr \^/ t \Q t ^Vabc{ma,mb,mc) {Aab,c-mc)n, 
Vl2a{ma)Sa{ma)V3ib{mb)Sb{mb) 7^^ ^ Vabcirric - Uc) 



Vabcirric) 



rirl 
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In words, this identity implies that after summing over n^ and Ub in (97), all the factors 



associated with nia and m^ become exactly what they should be according to the diagram- 
matic rules. Furthermore, the factors associated with nic are 77 — -, — 7 times exactly what we 
would have started with if we had considered the diagram with external lines 1, 2, 3 and 4 
stripped off. We have already seen in the previous section that a five-point function satisfies 
the rules; repeating that analysis here, the factor VahAjnf) of the five-point diagram cancels 
the 



1 



Vabc(mc) 



in the identity above to give exactly the correct result. 



If we had started with an even larger diagram, we would just apply eq. (98) repeatedly, 



reducing at each stage to the needed factors times a reduced diagram. To be a bit more 
explicit, consider a general diagram M„, factorized along a 'central' bulk to bulk propagator 
with a vertex with dimension Ae, taking the form 



M„ 



E 



S'e(me) 



VTif. 



EE n 

{nj {rrii} Vprop 



Si{mi -ni)\ {I^yz.e + my + m^-m. 



e/nie—ny 



6i - nii 



{me -riy -n^)! 



(99) 



n^'.: 



■.vert 



k{mi -ni,mj -nj,mk - Uk) 



{Aij^k + mi + TTij - rrik] 



n^—ni—rij 



{Uk 



rii 



rii 



X Rr^ 



where y, z are the lines connected to e on the left. The first product inside the brackets is over 



propagators, and the second is over all vertices. The identity (98) now says that, working 



"outwards in" according to figure [8} the sums on the n^'s exactly turn the Pochhammer 
symbols, shifted vertex factors, and shifted propagator factors associated with rrii into the 
appropriate final factors, times factors that allow us to consider a reduced diagram with the 



i legs stripped off. Thus, the poles of the factorization formula (68) are exactly given by the 
diagrammatic rules we have presented. 



A crucial aspect of the identity in equation (98) is that while it may seem complex, it 



is localized to a particular vertex in the diagram. It says that if we work from the outside 
of a tree diagram inward, toward the central propagator with an associated 6lr, then the 
sums over the rii variables simplify at each vertex independently. Without this property, it is 
hard to imagine that these diagrammatic rules could consistently reproduce our factorization 
formula. 

The relative simplicity of this localization allows us to immediately write down and check 
the necessary identities for any scalar theory. In a theory with both 3-pt 4'i(f)j4>k couplings 
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Figure 8: This is the general diagrammatic identity which is necessary to prove that our factoriza- 
tion formula and our diagrammatic rules are identical. The Pochhammer factors localize at vertices 
in the diagram, and then the sum over mj from the factorization formula can be performed locally 
on the diagram, working from the external legs inward towards the factorization or ^5lr propa- 
gator. Applying this identity to the sums turns the factorization formula into the diagrammatic 



rules. The precise form of the Pochhammer symbols is schematic, refer to equation (100) for the 
full form. 
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and an {R + l)-pt coupling, we would find a generalized identity 



i=l na,c,=0 



x14i3,...,aK3x(ma,3 -na,,]mx -nx) 7 x — — ^^^ (100) 






i=l 
.. Ki3,...,afl3x(0; nix - ^x) (Aai3...aa3X - mx)nx 

Vai3,...,an3x{0;rnx) nx! 

The first line of this identity is a product of all the vertices, propagators, and associated 
Pochhammers that lead into our {R + l)-pt vertex, the second line is the {R + l)-pt vertex 
and its associated Pochhammer. The third line is a product of simplified vertices with the 
n dependence eliminated, so that they take precisely the form that the diagrammatic rules 
dictate. In the final line we have an 'on-shell' rrij = version of the [R + l)-pt vertex, 
which will feed in naturally to the next vertex, leading towards the factorization propagator 
Slr- We have checked numerically that this identity holds in the case of theories with 3-pt 



and 4-pt interactions, and we have written out the general 4-pt vertex in appendix C.2 In 
theories with many different interaction vertices, there are equivalent identities involving all 
combinations of the vertices. 

5.2 The Diagrammatic Rules Satisfy the Functional Equation 

Next, we will show that the diagrammatic rules we have just presented satisfy the functional 
equation. As we discussed in section 2, the general functional equation takes the form 

Mo = {6LR-A)id-A-6LR)M+ Y, {^a^SbjM - 6ajSMM:^;l;; + 6a,S,,M^;i^^ (101) 

ab<k<ij 

where we recall that the shifted amplitudes are M^- '^j = M{Sab + 1, Sij + 1, 6ai — 1, Sbj — 1, ...), 
and we have chosen a specific propagator 6lji, which divides the Mellin amplitude into an 
M^ and M^_^^j^^. A feature of our diagrammatic rules is that they only explicitly depend 
on the left and right 5ij, and not those with i < k < j. Now the shifts such as M°/l* in 
the functional equation only involve L-R combinations, so when we plug in our factorization 
formula, we find that M^^'^* = M. This means that the first and second terms in the sum in 
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equation 35 cancel, so the functional equation reduces to 



Mo = {SLR-A){d-A-6LR)M+ J2 ^-bSijM^;^. (102) 

ab<k<ij 

Now, since the large Sij behavior of the diagrammatic rules are easy to read off, we can 
immediately see that the term linear in Slr in the RHS of this equation cancels. The 
constant term is by definition independent of 6lr, so this will have to be shown to match 
Mq on the LHS. Our strategy will be to first show that the RHS is independent of 6lr- But, 
any possible dependence on 6lr in the RHS is clearly in the part that falls like 6~j^. So, if we 
just want to prove that the RHS is independent of Slu, then all expressions at any step can 
be evaluated on their poles. Afterwards, we will show that the remaining (5LR-independent 
piece matches Mq]^ 

First, note the effect of the shifts in M^-'^L By assumption, ab and ij are on opposite 
sides of the 6lr propagator, so Slr — > Slr — 2. The other propagators are shifted as follows. 
Consider a propagator 6i on the left part of the diagram. If a and b are both to the left 
of this propagator, then Si — t- (5j — 1 in M^-'^^-, otherwise it has no change. A symmetric 
statement holds for propagators on the right part of the diagram. 

At this point, we can make effective use of the results from the previous subsection. 
Specifically, let us write the Mellin amplitude as 

M = y2MLi6ab,m,)^^^^MniSij,m,), (103) 



and rewrite the second term on the RHS of (|102|) as 

(5ab)n., 



b<k<ij me ^ ^ \E"ab=l 



nab,me) ■ ■ ' ■ "' 



ab<k<ij rUe '^ '^ \y^ran()=l a<b<k 




y^E'T'»i=l k<i<j 

where we have re-written the Sab^ij factors as Pochhammer symbols. Now, we can further 
simplify the terms in brackets using the result from the previous subsection that 



MLiSab + nab, me) = J^ | JJ afe + ^^^"1. j pj^^^^^ ^ ^^^ ^ ^>J^ 

n\=me \a<b<k '^^' / 



(105) 



^ We will not prove here that the overall numerical coefficient matches that in Afg, although we will 
obtain the correct parametric dependence on the couplings. 
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when all 6abS are evaluated on the poles. Substituting this into eq. (104]^, we find that 

poles ™'= 



ab<k<ij 



ab,ij 



(5p. -m^-1 



MiiSab + nab, rrie + l) JJ 



(Sab) 



nab 



a<b<k 



Uab^- 



X 



k<i<j 

rrieime + Ag - h)M\ 



iS: 



I] )rn. 



n- 



IJ- 



poles 



I poles ' 



(107) 



where in going to the last line we have shifted the summation index rrif, — )■ rrif, — 1 and used 



the fact that 



Se.{me-l) 



Ae+me—h 



Now, since {6lr- Ae){d- A^- Slr) = -4:6e{6e + A^- h), 



the poles in the above expression clearly cancel those in the first term in eq. (102)! 

We still need to show that the left-over ^L/j-independent piece, call it Mred, actually has 
the correct dependence on the remaining 5jj's to match Mq. To show that we produce the 
correct Mq, we can simply substitute Mj-ed into the functional equation for any of the other 
legs, and keep doing this until we have acted on every propagator in the diagram. After the 
last propagator is removed this way, the resulting Mq will be that for a contact interaction, 
which is simply independent of all the 6ij's. Thus, if we can prove that this procedure reduces 
the Mellin amplitude constructed from the diagrammatic rules to just a constant, then we 
will be done|3 

However, this follows immediately if we use a nice property of the relevant conformal 
Casimirs. The RHS of the functional equation (35) is the action of the Casimir (YIi^^l JiYi 
where L denotes the set of all external lines to the left (with respect to an a priori chosen 
vertex at the far right of the diagram) of the propagator being acted upon. Because of the 
tree structure of the diagrams, for any two propagators with corresponding L, L', either L 
and V will be disjoint or one of them will be a subset of another. In the former case, it is 
obvious that the two Casimirs commute. A short computation shows that they commute in 
the latter case as well: 

2 / \ 2' 



E-'. AT, J' 



jeL / XieL'cL 
^It is necessary here to use the fact that 






E^. ■ E-'. 



ieL' 



Kiel' 



0. (108) 



E E 



n 



rUe a<b<k 



nab'- 



E 

a<b<kJ2^ 



E 



y-j- [Sab)n'^^ „ 



-rrie+l a<b<k 



n 



ab 



(me + 1) ^ n 

Y, n^'^=mc + l a<b<k 



b 
iSab)n'^ 

n" 



(106) 



^ We will address quite generally the issue of homogeneous solutions in appendix A 
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Therefore, all the Casimirs associated with propagators commute with each other. Let us 
act on M with the Casimirs corresponding to all the propagators in the diagram. Since the 
Casimirs commute, we can act with any one of them first, so the result must be independent 
of all the (5jj's, which completes the proof. 

6 The Flat Space Limit of AdS/CFT and the S-Matrix 

One can obtain a holographic description of the flat space S-Matrix [271 EHl EHl EQl El [3ll |32l 
155] by taking a limit of the AdS/CFT correspondence in which the curvature of AdS goes 
to zero. This fact has received more attention recently, and in particular in [3] one of us 
argued for an intimate connection between the Mellin representation of AdS/CFT correlation 
functions and the fiat space S-Matrix. 

Roughly speaking, the claim is that if one interprets the 5ij as the kinematic invariants 
Pi-Pj of a scattering amplitude, then in the large 5ij limit the Mellin amplitude will reproduce 
the fiat space S-Matrix. This result has a very simple physical justification. Time translations 
in global AdS are generated by the dilatation operator of the CFT, so the dimension of an 
operator in the CFT directly translates into the energy of a corresponding bulk state. Since 
the AdS radius -R is a dimensionful quantity, what we really mean by the fiat space limit is 
ER — )■ oo for all bulk energies E. We have seen repeatedly that the Mellin space coordinates 
5ij are related to the dimensions of operators in the CFT, so it is natural to expect that the 
large 5ij limit of a Mellin amplitude is related to the physics of bulk states with large ER, 
which compute the fiat space S-Matrix. We will now show how this works quantitatively. 

6.1 The Flat Space Limit of AdS/CFT Factorization 

At a computational and pictorial level, we know that Witten diagrams describe scattering 
events in AdS spacetime. In [3j, an explicit formula was given relating the large 5ij behavior 
of a Mellin amplitude to the corresponding fiat space scattering amplitude. With the present 
normalizations, this formula reads 

h " n roo 

M{Sij) ^ y nf(^ /, ci/3/3^S^-^-ie-'^T(p, -p, = 2/3(5,,) , 6,, » 1 . (109) 

where the fiat space S-Matrix is 5* = 1 -|- i{2'7iY~^^6 (J^Pi) ^• 



The transformation in (109) may appear complicated but its effect is relatively simple. To 



understand why something like ( 109 ) is necessary, first consider the case of massless scalars 
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scattering through a single contact interaction with some fixed number of derivatives. In this 
case, the only dimensionful parameter is the coupling constant, which is an overall factor 



in the amplitude T{pi -pj). Thus, (3 in (109) just factors out by dimensional analysis, and 
the Mellin amplitude is simply proportional to the fiat-space scattering amplitude. The only 
effect of the transformation is to modify the overall coefficient. The reason that a tran- 
formation is needed for more general theories is that this overall coefficient depends on the 



dimension of the interaction. It was shown in [3] that for all scalar contact interactions, (109) 
produces exactly the correct overall coefficient in order to match the Mellin amplitude. This 
is a very general check, since one may typically consider an arbitrary scattering amplitude 



as a linear combination of such interactions below some cut-off. So, roughly what (109) says 
is that to get the flat-space S-matrix from M{6ij), one simply performs a series expansion 
at large 6ij and goes through the series term by term, altering the coefficients by hand in a 
way that depends only on the power of the 6ij's. 

Now, the goal of this section is to show that the fiat space limit of the factorization 
formula 



M{6,,) 



y. 47r^r(A -h + l)m! L^(^,,)i?^(,5,, 

m=0 



{A-h + l)m 5lr - A 
reduces to the usual factorization of scattering amplitudes. 
We recall that Lm is given by 



2m 



(110) 



I^m\Oi 



my-'i] ) 



Y, M^(5., + n,,)n 

J2nij=m 



(S. 



IjJUi 



Kj 



n. 



ly 



(111) 



=A+2m 



and similarly for Rm- We start by considering the limit of 5ij ^ 1 with fixed internal and 
external dimensions A and Aj. In this case, the fiat space limit will give rise to scattering 
amplitudes of massless particles because the mass squared of the bulk fields is of the order 



of the AdS curvature. This turns equation (110) into 

47r'^r(A-/i + l 



M{5,, 



ijj 



>LR 



E 

m=0 



771! 



(A - /i + 1), 



-Lm\Pij)H„i\pi-j 



^J) 



(112) 



We see that the multiple poles of the Mellin amplitude at Sir = A + 2m all contribute 
to the pole of the scattering amplitude at (Xli-i "PiY = 0. In order to show that the sum 



Ji=l 



over all contributions reproduces the correct residue of the scattering amplitude, we need to 
understand the large 5ij limit of L^. It turns out that at large (5j,-, L^ simplifies significantly. 



Lm\0-: 



m\^i] ) 



1 



^l(j:^.-^)-iML{t5, 



5ij > 1 



(113) 



t=i 
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We will now prove this relation. To that end, consider the example 



^^(5..) = n(^«^-)'^«. 



(114) 



i<j 



This set of functions is very broad and can be used as a basis. Thus, to prove (113) in 



general, it will be sufficient to prove it for (114). This is a convenient basis to use because 



for this type of left Mellin amplitude, one can perform the sum in Lmi^ij) explicitly. 



ETT ("»i + (^ij)nij ,^ . 
11 „..l y'^ijJH 



(E ^ij + E «' 



SLR=A+2m 



i] )m 



m\ 






ij Jaij 



SLR=A+2m 



(i(EA,-A) + E«.,-m), 



m! 



Kj 



I] ) aij 



(115) 



where, in the second line, we used the identity (|80|) and, in the last line, we evaluated 

k 
i<j 



Ef<j ^ij = (ELi ^» ~ '^Lij)/2 on the pole 6lr = A + 2m. For large 6ij, we have 



M\6,,) ^ l[{S.,r^ , LUk 



Q(EA.-A) + E 



an — m) 

'■J ' m 



Kj 



m\ 



n(^^^)"''' (116) 



Kj 



in perfect agreement with (113) 



Let us now substitute this simplified form into (112) and invoke the fiat space limit 



formula (109) for M^ and M^. We obtain the large 5ij limit of the Mellin amplitude in 



terms of the left and right scattering amplitudes: 



M{5, 



vr'^ A C; 



4(5 



\n n r(A,) 7o 



dh 4^^^^'^' ^'^-'- V^(p, ■ p, = 2hk^) 



c^/3i?4^^^^^>'^^'^'''V^(p, -p, = 2/3^5,,) 



;il7) 





oo 



E 

m=0 



d d 
m\ \dtLdtpi 



{tLtRf-^-^e-P^/'^-P^''^' 



r(A-/i + l + m) 



tL=tR = l 



where we have rescaled the integration variables /S^, Pr. This looks like a rather complicated 
expression but it simplifies dramatically due to the following identity 



E 

m=0 



1 ( d d Y {tLtRf-'^-^e-^^"^-^^"^ 



ml \dtL dtji 



T{A-h + l + m) 



(3t^e-^'Si(3L-(3n) 



118) 



tL=tR = l 
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which we prove in Appendix C.3 Finally, the large Stj behaviour of the Mellin amplitude 
simplifies to 



M{5,, 



ij) 



h " 

2 Alr(A, 



TT- -r-r Ca^ 



d(5[5 



lyA^-h-i-p T^JPi-Pj = 2l36ij)T^{pi ■ pj = 2l36i. 



H) Jo —2(36lr 

in agreement with the factorization of scattering amplitudes for massless scalars, 

T^{pi-Pj)T^{Pi-Pj) 



(119) 



T{Pi ■ Pj) 



iEi=lP^ 



(120) 



6.1.1 Massive Propagators in Flat Space 

The mass of a scalar field in AdSf^+i is A(A — d)/R'^, where A is the conformal dimension 
of the dual operator and R is the AdS radius of curvature. Thus, in order to keep a finite 
mass in the fiat space limit i? — ?■ oo, one must scale the dimensions A — )■ oo of the scalar 
operators. 

More precisely, one considers the limit of large A and Sij with fixed ratio Sij/A"^. Let us 
then return to (110) with (113) and (109) applied to M^ and M^, and study this limit 



M(5i 



IT 



n 



Lr(A 



ij JO 



I 1 



rf/3ij/3|^«^' ' V(pi-p,=2/3K5.,) 





oo 



E 

m=0 
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In appendix C.3 , we prove that the limit of the last line is 

2/3l 



SWl - PR)e-^'P'l 



2I3l5_ 



LR 
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'1221 



Thus, the limit of (|121| is 
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dl3(3-' 



-.T.^,-h^i.~p T\p,-p, = 2P6,j)T^{p,-p, = 2^6.. 



2 A^A r(A,) 7o -2PhR + A2 

as expected from the factorization of scattering amplitudes 

T^{pi-Pj)T^{Pi-Pj) 
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nP^■PJ) 



{EliPi)' + M' 



(124) 
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where we identified A^ with the mass squared of the exchanged particle. 

It is worth noting that the poles of the MeUin amphtude at Sir = A + 2m turned into 



the pole at 5lr = A^/(2/3) in the integrand of (123). The reason is that, in the flat space 
limit, the infinite sum over poles is dominated by the poles with m of order A^ ^ 1. More 
precisely, in the limit of large 5ij with fixed bijj Is^ ^ the infinite sequence of poles of the Mellin 
amplitude gives rise to a branch cut along the positive real axis of bi^RJ t^ . 

6.2 Diagrammatic Rules 

In the last section, we have shown that the flat space limit of the AdS factorization formula 
reduces to the usual factorization of scattering amplitudes. In this section, we will show that 
the fiat space limit of the AdS Feynman rules proposed in section |5] gives the usual Feynman 
rules for scattering amplitudes. 

The fiat space limit corresponds to the large bij behaviour of Mellin amplitudes. If we 



take this limit with fixed A's, then equation (92) simplifies to 






r 1 1 ^ 



where the index a = 1, . . . ,s labels the internal propagators of the Witten diagram. The 
function M(mi, . . . ,ms) is computed using the diagrammatic rules of section [s] This large 
6ij behaviour of the Mellin amplitudes should be compared with the prediction from the 



flat space limit formula (109). Inserting the appropriate scattering amplitude for massless 



scalars in (109), one obtains 




2^ ' '■':iSa 

1=1 / a=l 

where we have set to 1 the coupling constants associated to each interaction vertex. We 
conclude that the dependence on the kinematic variables 6ij is the correct one. To finish the 
proof we just need to show that the overall normalization also agrees, i.e. we must show that 

J2Mi„r„...,,n.)^^-l(f[^)r(l±A.-H-^^^. (127) 

{ma} \i=l ^ *V \ i=l / 

We will perform the sum over {nia} recursively, starting from the bulk propagators closer 
to the external legs of the tree level Witten diagram. 
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The first case to consider is a part of tlie Witten diagram tliat connects one bulk propa- 
gator (Ag, mg) to 3 external legs (Ai, A2, A3) like in the left diagram in Fig. [s] In this case, 
we need to compute the sum 



J^ VAiA2A7(0,0,m7)S'A7(m7)VA3A7A8(0, 



mj,ms) 



(128) 



m,7=0 



_ -27r^Ai27A378r^(A7 -h + l)T ( Ai+a^+As+a^ _/,_!) 

"m8!(Ai23,8 - l)-msr (^l+f±^ -h)T (A.i±f±A. _ h) ' 

It is then easy to use this result twice to compute the sum 

00 

y^ ^AiAaArlO, 0, 7717)8 Ar{rnr)VA,^ArAs{0, "^7, m8)SAs{ms)VAiAsAg{0, mg, ITlg) 
r?i7,mg=0 
_47r2'^Ai27A378A489r2(A7 -h+ l)T'{As -h+l)T ( A1+A.+A3+A4+A, _ ^ _ 2) 



(129) 



m9!(Ai234,9 - 2)_„,r (Al+A^+Az _ /,) r (Aa±Az±A, _ ^) ^ ( A4+AB+A. _ ^) ' 

corresponding to the 2 leftmost bulk propagators in the left diagram in Fig. |5] By using this 



rule recursively, we can compute the sum in (127) for some Witten diagrams. However, the 



general Witten diagram requires another type of basic sum, corresponding to a part of the 
Witten diagram connecting 4 external legs (Ai, A2, A3, A4) to a bulk propagator (A9,m9) 
like in the right diagram in Fig. |5| 

00 

^ ^AiA2A7(0,0,m7)S'A7(m7)VA3A4A8(0,0,mg)S'A8(mg)VA7A8A9('^7,'^8,"^9) (130) 
mr ,ms=0 
_ 47r^^Ai27A348A789r^(A7 -h+ l)r^(Ag -h+l)T ( A.+A.+A3+A,+A, _ ^ _ ^^ 
m!(Ai234,9 - 2)-^r (A1+42+A7 _ ;,) r ( A3+A4+A8 _ ;,) r (Ai+As+Aa _ /,) " 



Using (128) and (130) recursively, it is easy to show that a general tree level n-point 



Witten diagram with s internal propagators has 
^M(mi,...,m,) = r('^^ 

{ma} ^ 

where 



ft-p)(-.''rnr'(A.-A+i)n-7^ 

' a=l vert J- I 
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ijk 
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A., = ?r(^i±4i±^-A)n^- 



r(Ai 



+A2+A3 



h) 



Lr(A,) 167r2'^ntir(A,-/i + l) 



131) 



is the 3-point Mellin amplitude. Finally, using the fact that a tree level n-point Witten 
diagram has s = n — 3 bulk propagators and n — 2 cubic vertices, it is straightforward to 
prove (127). 
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7 Discussion 

We have argued that Conformal Field Theory correlation functions have a natural home in 
Mellin space, and we have given dramatic evidence for this claim in the case of CFTs with 
a weakly coupled AdS dual. 

The easiest way to summarize our results is to list the profound similarities between Mellin 
space for AdS/CFT and momentum space for scattering amplitudes. As with scattering 
amplitudes in momentum space, in Mellin space 

• CFT correlation functions have poles corresponding to the exchange of operators, which 
we have dubbed 'OPE factorization', and on these poles the correlation functions 
factorize into lower point correlators; 

• the differential equations that define AdS/CFT correlators as Green's functions turn 
into simple, purely algebraic functional equations for the Mellin amplitude; 

• there are simple diagrammatic rules that enable a direct construction of the Mellin 
amplitude corresponding to any Witten diagram; 

• the Mellin space 'momentum' flows through these Witten diagrams in such a way that 
it is conserved at all vertices. 

Furthermore, the connection becomes totally explicit when we take the fiat space limit of 
AdS/CFT, where Mellin variables turn into flat space kinematic invariants via 6ij — )■ Pi ■ pj 
and the Mellin amplitudes themselves reduce to the fiat space S-Matrix of the bulk theory. 

In this paper we have only dealt in a precise and systematic way with CFTs dual to 
theories of scalar fields living in the bulk of AdS, and we have only computed the correlators 
at tree-level, i.e. at leading order in 1/A^. We expect that this is only the beginning. It will 
be interesting to understand how our results generalize to theories with higher spin fields, 
at loop level, and to more general CFTs that do not have a weakly coupled AdS dual. Let 
us conclude with a few comments about these possibilities. 

As discussed in appendix [B| there exists a straightforward method for generalizing our 
results to vector and tensor fields. However, more efficient methods might very well exist 
which more naturally incorporate the helicity structure of the fields dual to the conserved 
currents and the stress-energy tensor. In this respect, (standard) momentum space seems to 
have an advantage since the current conservation condition, which is a differential equation 
in position space, turns into the more tractable algebraic equation p^J^ = in momentum 
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space. Indeed, recently some simple results have been obtained for tensor correlators in 
momentum space [H] . Second, computing gravitational amplitudes using Witten diagrams is 
tedious because gravity has an infinite number of interaction vertices. To efficiently compute 
these amplitudes in Mellin space, we would require some version of the BCFW recursion 
relations [IHl HHl HZ] analogous to the one developed in momentum space for AdS |1H1 SH] • 
Unfortunately, conformal invariance is obscured in momentum space. Moreover, in some 
cases (logarithmic) divergences in the Fourier transform of CFT correlation functions in 
momentum space lead to conformal anomalies which complicate the analysis for higher-point 
functions, see for example [50]. These issues do not arise in Mellin space. 

At loop level in the bulk of AdS, there may be two distinct but natural forms for the Mellin 
amplitude to take - one generalizing our factorization formula, and the other generalizing 



our diagrammatic rules. As we discussed in section 2.1 , we expect that at loop-level and even 
non-perturbatively, the Mellin amplitude will be meromorphic, with poles corresponding to 
the exchange of operators in the OPE. So at loop-level it should still be possible to write 
the amplitude in a form that makes the poles corresponding to various multitrace operators 
manifest; a form with exactly these properties was found for the 4-pt amplitude in [3j. 
However, given the success of our diagrammatic rules, we also expect to find loop-level rules, 
where infinite sums replace the momentum space integrals familiar from fiat spacetime. In 
these diagrams, factorization may be obscured, but perhaps bulk unitarity will be manifest. 
It is not hard to guess a generalization of our rules from tree to loop level, but we will leave 
the exploration of their nature and their physical interpretation to future work. 

We should note that the complete holographic computation of correlation functions in for 
example N' = 4 SYM using supergravity is more involved than just the evaluation of the Wit- 
ten diagrams themselves. In particular, in order to obtain for example the five-dimensional 
Lagrangian a Kaluza-Klein reduction has to be performed [51] which beyond the first order 
is rendered computationally difficult because of the non-linear relations between the Kaluza- 
Klein coefficients and the eventual five- dimensional bulk fields [52]. Using superconformal 
invariance and inspired by the results obtained in this way, a general conjecture for certain 
scalar four-point functions of A/" = 4 SYM in the supergravity limit was written down in 



Results like those of [21] for AdS/CFT four-point functions were often conveniently ex- 
pressed in terms of the so-called D-functions which were first introduced in [16]. Such 
D-functions correspond to contact Witten diagram in AdS and as we have seen in equation 



(39) their Mellin transforms are just constants. One may therefore attempt to use existing 
techniques in the literature [16j to convert bulk to bulk propagators to contact interactions 
and amplitudes to sums over D-functions, and then reformulate these results as Mellin am- 
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plitudes. We have verified that this leads to answers in a different form than we presented 
above and the two forms are related by resummation (as, for example, in the four point 
calculation of [3]). The results of this paper indicate that the Mellin transform may be a 
more natural object to describe higher-point correlation functions. We hope that this can 
enable an extension of the current results for e.g. A/" = 4 correlation functions in the su- 
pergravity limit to more than four external points. We also expect Mellin amplitudes to be 
useful in A/" = 4 SYM beyond the supergravity limit. An important question is if higher 
point functions of single-trace operators can be constructed from the knowledge of 2-pt and 
3-pt functions of single-trace operators. The factorization of Mellin amplitudes provides a 
concrete way to do this in CFTs whose AdS dual can be constructed from cubic interaction 
vertices. More generally, we can hope that a BCFW-like construction is possible for type 
IIB string theory in AdSs x S^, which recursively reduces n-pt to 3-pt Mellin amplitudes. 

Will our results extend usefully to general CFTs? The present definition of the Mellin 
amplitude [1] includes factors of T{6ij) which are especially convenient for theories with a 
perturbative expansion, where the existence of operators with dimensions Ai and A2 implies 
the existence of an operator with dimension A1+A2, up to perturbative corrections. We have 
only studied theories with a small number of low- dimension operators, so it may be interesting 
to investigate Mellin amplitudes in say weakly coupled A/" = 4 SYM (not necessarily at large 
A^) , where the perturbation expansion involves a larger number of operators. But the larger 
question is whether the Mellin representation can be useful at a non-perturbative level. Our 



arguments based on the OPE from section 2.1 suggest that it will be, but we will need 



to make those arguments more precise in order to apply them quantitatively. It would be 
interesting to investigate this question by studying simple, exactly solvable models such as 
minimal models. One by-product of our analysis that could be immediately useful is the 
functional equation, since its homogeneous solutions are nothing other than the conformal 
blocks. In the future we expect to present results using the functional equation to obtain 
explicit expressions for the conformal blocks. 

Finally, we have extended the discussion of [3] , showing that when we take the flat space 
limit of a classical scalar fleld theory in AdS, the Mellin amplitude of the dual CFT morphs 
into the S-Matrix of the bulk theory. This suggest that Mellin amplitudes can be used to 
give a holographic and non-perturbative deflnition of the gravitational S-matrix, through the 
flat space limit of AdS. The main open question is how much progress can be made in the 
computation of Mellin amplitudes beyond the planar limit (i.e. tree level in AdS). It will be 
exciting to pursue this research avenue towards a holographic description of flat spacetime. 
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A The Functional Equation and Large Sij Behavior 

We showed in section[5]that the diagrammatic rules are identical to the factorization formula. 



and also that the diagrammatic rules satisfy the functional equation of 2.3 Furthermore, we 
know that our formulas reproduce the poles of Witten diagrams from the discussion of the 
factorization formula in section |3} so if there is any difference between our formulas and the 
Witten diagrams then this difference must be an analytic function. However, by linearity of 
the functional equation we find that any possible difference should also satisfy the general 



functional equation (35) with Mq = 



Our formulas explicitly vanish as 6ij — )■ oo, and on physical grounds we expect that 
Witten diagrams must be polynomially bounded in this limit. Among other things, if this 
were not true then our Mellin amplitudes would not reproduce the flat space S-Matrix, as 
we know they do. So to really prove that our formulas are identical to the Mellin amplitude 
of the corresponding Witten diagrams, it suffices to prove that no polynomial can satisfy the 
homogeneous functional equation. 

Let us see why a pure polynomial cannot satisfy the homogeneous functional equation. 
Roughly speaking, one would expect this to follow because the functional equation requires 
a certain periodicity, and polynomials clearly cannot be periodic, but we can easily make a 
more precise and direct argument. The general homogeneous functional equation is 

i5LR - A) (d - A - 5lr)M + J2 {Sa^5b,M - 5a,5uM'^^ + SabS^.M^f^ = (132) 

ab<.k<ij 
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Let us use a basis for the 5ij that is unconstrained by identities involving the dimensions Aj. 
Any polynomial in this restricted set of 5ij variables will have a term of highest degree D of 
the form 

n 



H{S^,)= Yl f(d^^)IlK 



(133) 



Y:d^j=D 



i<j 



Plugging this into the functional equation, we see that the term in that equation of greatest 
degree is unaffected by the shifts of the 6ij, so the first two terms in the summand cancel 
and we find that 



i5LR-A){d-A-SLR)+ Yl ^^b^'- 

ab<k<ij 



H{6, 







(134) 



The term in brackets of order (5|^ will automatically cancel, but the term of order 6lr will 
not, at least for generic (physical) values of the Aj. It also cannot be canceled by the next- 
to-leading degree terms in M because the order 6fj coefficient of those terms also cancels, so 
the only solution to this equation is H = 0, completing the proof. 

Now, for completeness, we will present an argument that the functional equation can be 
directly used to argue that our factorization formula correctly computes Witten diagrams. 
(Notice that in the body of the text we proved this via a detour, through the equivalence 
of factorization and the diagrammatic rules.) Namely, we will prove that our factorization 
formula satisfies the functional equation when it is evaluated on its poles in 5^/?. We know 
that the poles and residues of the factorization formula are identical to the poles and residues 
of the Witten diagram in question, from the development of section [3} and so from the 
argument above we then obtain that our factorization formula is indeed correct. 

We want to prove that our factorization formula is sufficient to determine the Mellin 
amplitude. More precisely, we claim that 

Res(m) 



M = E 



m=0 



6lr- A- 2m 
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Resim) 
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Rm (136) 



ALR=A+2m 



where of course R^ takes the same form as the left piece, so we have not written it out 
explicitly. The arguments from section |5] immediately apply, so the functional equation 
reduces to 

Mo = {5LR-A){d-A-5LR)M+ Y Sa^SijM^^g (137) 



ab<k<ij 
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The expression on the right hand side appears to be of order 5\p^M as 5lr — )■ oo, but as 
before one can check exphcitly that this dependence cancels between the first term and the 
summand. Then the leading behavior for large 5lr is of order 6lrM, which approaches a 
constant, independent of Slji. Thus up to a constant, we can evaluate the right hand side of 
this equation on its poles. We will find a dramatic simplification. Let us focus on a particular 
term in the summand and study 

AB ij Ai,Bj Z^ SLR-A-2m-2 T{A-h + l + m) ^-^ '"^ ^^ ^-^ ' ^ ' 



E in*' 
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ijjriij 






{Sab + i; 



nAB 
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nAB^- 



SABM^{Sij + riij] Sab + riAB + 1) 



where we have capitalized the AB, IJ to differentiate them from the ij that are being summed 
over, and again we are leaving right piece implicit because we will be manipulating it and 
the left piece in an identical way. Furthermore, since we are ignoring the constant piece and 
only consider the the poles and their residues, we are implicitly evaluating all the Sij in this 
expression at the pole Slr = A + 2m, in particular the Sab and Su are now also assumed 
to be subject to this constraint. Focusing on the second line, notice that 

^-^^^^^^Sab = {UAB + 1) f^^^r^l (139) 

uab^- {nAB + ly- 

Let us use this fact and switch the order of summation, so that we sum over the AB labels 
inside the sum over the rijj, 

J2 E (UAB + 1) n ^ ^^^%^M^(5., + n,,; Sab + uab + 1) (140) 



'ynij=m A<B 






\ 



The result is that we obtain precisely the (m + l)th term of the series, except for {hab + 1) 
type factors. However, when viewed from the perspective of the [m + l)th term, this just 
means that we must multiply by 'Yl,^ij = rn + 1. This follows from the following identity 
regarding sums over partitions of integers, 

k 

Y^ E(n^ + l)F(ni,...,nA + l,...,nfc) = (m + l) ^ F(ni,...,nfc). (141) 



Thus we find that the entire summand in (137) has the following poles 



J:^rn{2^-d + ^-^) ,^!'lT\^ (142) 

m=0 
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Evaluating also the first term on its poles gives 



Res{m) 



{5lr -A){d-A- 6lr)M ^-T 2m(2A -d + 2m) ^^ , (143) 

^0 dLR-A- 2m 

proving that the right hand side of the functional equation is analytic in 6lr- To complete 
the proof that our ansatz solves the functional equation it is sufficient to show that the right 



hand side of (137) tends to Mq as Slr -^ oo. This depends on the theory and the structure 



of the diagram, and proving this was precisely the subject of section [5] 

B Exchange of Vector Fields 

In this appendix, we wish to examine how Mellin amplitudes factorize when a gauge boson 
is exchanged in a bulk to bulk propagator. This analysis can be generalized to the exchange 
of a higher spin field in the bulk. We will only consider amplitudes with external scalars. 

We need a formula that relates the bulk to bulk propagator to the bulk to boundary 
propagator 

G^BBA^^^y) = j X{c)dc / rfSG'^^_(^_^^)(x,2;)G'^'^(/,_^)r7p^, (144) 

where Gqb,a indicates the bulk to boundary propagator for a spin-1 field of dimension A, 
and Gbb,a is the bulk to bulk propagator for the same field. The existence of such a formula 
follows from the existence of the analogous formula for scalars. All we will assume about the 
function x(c) here is that it has a pole at c = A — h. (A formula of the kind that we need 
was developed in [53] but we will not need its detailed form.) 

By inserting this formula in a Witten diagram we get 

A{xi, ...Xn)= / X{c)dc / d'^zAl{xu ■■■Xm, z)A''j^{z, Xra+1, ■ ■ ■ Xn)r]^,y. (145) 

Let us lift the vector fields A^ and A^, to vector fields on the boundary of AdS in the 
embedding space. As explained in [SU ST], a spin-1 primary operator A^{x) is uplifted to a 
transverse spin-1 field Am{P) on the light-cone of M'^+^, such that 

apM 
P'' Am{P) = , A,{x) = -g^^MiP) , (146) 

where P^ = (1, x"^, x™) is the Poincare section of the lightcone. 

It follows that 

r^^^A^A^ = VmnAI'A^ , (147) 
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which uphfts equation (145) to the embedding space. 

Conformal invariance tells us that A^ must have the form 

m 

A^{P„ ...P^,P) = Y, ^liPi, ---Pm, P)P^ (148) 

p=0 

with Pq = P and where the A^j^{Pi, . . . P^, P) are m + 1 scalar functions. 

We will also need the fact that the function A^ has a conformal weight at the point Pq 
which is given by 

^-^ \Aq + l,iiq = p ^-"^^^ 

where Ag are the dimensions of the external operators in the original amplitude A. 



Inserting this into (|145|), we get 
A[Xi, . . . Xfi) 



/p HI, II, 

X{c)dc d'^PY, Yl Al{P^,...P^,P)A%{P,P^+^,...Pn){Pp-Pq) 
•^ p=l q=m+l 

m n 

= E E ^P,(A,...,Pn)(Pp-P,). (150) 

p=l q=m+l 

Note that the sum over p starts from 1, because the terms involving P ■ Pi etc. have dropped 



out because of the transversality constraint (148). From now, the range of the sums will 
be kept implicit. Also in the term {Pp ■ Pq), one factor comes from the left, and the other 
factor comes from the right. In general, the effect on the Mellin amplitude of multiplying the 
correlation function by the factor Pi ■ Pj is to send M{6ij) — )■ 6ijM{6ij — 1). We can therefore 
for the moment consider the Mellin amplitude in the absence of this factor, and reintroduce 
it later. But, this is exactly of the form we encountered in scalar theories, where we have 
seen how factorization works. So, repeating the analysis for scalars, we now find that 



Mpq{6ij) = / xic)dcLj, X Rq, (151) 



where Mpq is the Mellin transform of Apq. Here, the subscripts on L and R come from (150) 
and do not indicate the number of particles inside L and R. We have kept that information 
implicit to avoid clutter. We have 

Lp = [[d6UdlU M,^(4, /.) n ^^^"^rl^" ~ ^"^ (152) 

where the constraints on the 6 are now 

5^4 = Af (153) 
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where A^ is given by ( 149 ) . 



We are almost done now, since we can repeat the analysis for poles when a scalar is 
exchanged. We find that Mpq contains a pole when 

J](<5,, + n,,) - ^(5^ Af - /. - c) = 0. (154) 

ij i 



However, since ^^ Af = ^^ Aj + 1, we can rewrite (154) as 

5LR = ^/\i-2'^5ij = h + c-l + 2m. (155) 

i ij 

In order to obtain M from the Mpg's, we take M = J2pq^pq^^pqi^pq ~ !)• However, since p 
and q are always on opposite sides of the propagator, their shift has no effect on the position 



of the poles in 5lr- Therefore, (155) is the position of the poles in M as well. If we now 
assume that the function x(c) has a pole at c = A — /i, we find that for the exchange of a 
gauge boson, the Mellin amplitude has a pole at 

Slr = A-1 + 2m, (156) 



with a residue that can be read off from (151) 



A very similar analysis can be performed for gravity. As we mentioned above, however, 
the primary complication in computing graviton amplitudes is that Witten diagrams involv- 
ing the interaction of gravitons are inordinately complicated. It would be interesting to see 
if this can be ameliorated using BCFW recursion relations and to compare their form to 
the recursion relations developed for correlation functions of stress tensors and conserved 
currents in 



C Some Technical Developments 

C.l Shadow Field Identities 



We saw in section 2.1 that to implement factorization and unitarity, for each operator with 
dimension A we need to introduce a shadow operator with dimension d — A. Our goal is 
to understand factorization in the Mellin representation, so we should first understand the 
relationship between the Mellin amplitude of a product of operators Oi{xi) and the identical 
amplitude where one operator, say Oi(xi), is replaced by its shadow Oi{xi). The result is 

M(4) = -fn^. f[dS]M{S.,) n ^^^-^-^^ffi - ^-^ . (157) 

^ ^' -^ l<i<j -*- \"i'j) 

53 



Now we will derive it. 

The starting point is an identity relating propagation with dimension A and d — A: 



{-2P-XY-^ 7r'^r(/i-A) 



2P ■ P'y-A (_2P' . X)^ 



1581 



where X is a point in the bulk of AdS, but P and P' are boundary points. If we apply this 
identity to the computation of an AdS amplitude A and then represent A using a Mellin 
amplitude M, we find 



A{xi, 



1 



I 'J^r, 



7r^r(/i- Ai) 



r(/i-Ai 



^yr^^^^^i^A,,...., 



[d6] 



(xi - y)2(<i- 



[d6]M{6,,) H 



I -^r, 






(159) 

n 



i<j 



Or in other words, we have derived the relation (157) between the amplitude M and the 
equivalent amplitude with Oi ^ Oi. 



Now, we would like to use this relation to evaluate the R piece in eq. (65). Let us rewrite 
it with variables relabeled in a way that suits our present purposes: 



Mfc+i((5, 



^3^ 



r(A, 



+1 



,, jmmim'-^^. 



(160) 



Z<J 



where the integration variables are constrained by X^i^^j ^«i ~ ^if ^^^ ^^^ amplitudes M{Sij) 
and M{6ij) differ only through A^+i = d — Ak+i in the first external leg (for all other legs, 
Aj = Aj). In general, this applies only in the case that the Sij variables also satisfy similar 
constraints, Yli=/=i ^ij ~ ^ii whereas in our factorization eq. (65) we need to be able to take 



5lr = A/c+i + 2m when m 7^ 0. So this identity is not yet directly usable. However, for the 
special case of three-point functions (A; = 2), the integrations are vacuous, and we can derive 
a shadow field identity for arbitrary 5ij. To do this, let us define new variables 5ij that do 
satisfy the constraints, and write 



M3(5, 



1 T{5i2 - ^12)r((5; 



12 



r(A3-/i)r(5i2-<5i2)r(5i2) 
r(5i2-Mr(4) 



T{5ij)T{5ij - 5ij) 



[d5]M{5,,) n 

i<3 r((^ij-, 



r((5i2-5i2)r(5i2) 



M3((5, 



(161) 



One can now substitute this relation into eq. (65), with (1, 2) — )■ (n — 1, n) for continuity of 
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index labels and A3 — )• c + /i, ^13 — > In-i, (^23 ^ ^n, and constraints that give 



2Sn-l,n = A„_i + An- h + C 

25n-l,n = A„_i + An- h- C 

2ln-l = h + C + An-1 - An 

2ln = h + C + An — A„_i 

25n-l,n = A„_i + An — Slr. 

Thus, we find that on c = A — h, 6lr = A + 2m, we have 



(162) 



R = -r(A - h) ^}^ ,^^^;^ M,^{6n.,,nJi) = -^^-^^^-^M,{Sn-l,n + m){Sn-l,n)r 



/ A„_i+A„-A 



{A-h+l] 



This is sufficient to prove the factorization formula (43) for all poles of 6lr in n-point 
functions when adding on three-point functions: 



(h-A). 



A„_i+A„-A 



-M^iSijJi) 



(163) 



Now, the key to generalizing beyond three-point functions is that when we factorize the 



propagator to obtain eq. (65), either side of the diagram may be chosen to be the "Left" or 



the "Right", and the answer must be the same regardless. Instead of adding a three-point 
function onto a fc + l-point function, we could have made the fc + 1-point function the "Right" 
piece. 



Applying eq. (65) in this way with L ^ R, we obtain 



Mn{Si,) ~ Yl 



-^^{h-A) r^ ^.miSn-l,n)r 



6lr- A- 2m 



Ms^{-iy 



X L. 



m! 



(164) 



A„_i+A„-A 



By the constraints, we have Sn-i.n+iTT' = " ^ 2 " — ^^ ^^^ above formula. Matching residues 

therefor 

n„ i^ij)n., \ (h - A)_, 



of the poles in 6lr in eqs. (164) and (163), we therefore obtain 

k 



Y^nij=m 

{h-A){-ir 



'.i<j 



riijl I I A„_i+A„-A 



m! 



A„_i+A„-A 



X L. 



(165) 
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After some simplifications, this reduces to 
T{A~h)m\{-iy 



{A-h + 1] 



J2nij=m 



Hi 



which proves the identity (66) 



n 



{hj) 



IJJUij 
77 . .1 



(166) 



C.2 General Diagrammatic Identities and Vertices from the Func- 
tional Equation 



We saw in section 5 that although the factorization formula (43) looks diagrammatically non- 



local, in that the residues of the poles in one 6i appear naively to depend on the positions 
of all the other 6/s, it in fact satisfies a set of diagrammatic rules for each vertex and 
propagator that depend only on adjacent pole positions. A key technical aspect of this fact 



was an identity (98) that allowed the rules to be proved recursively, so that larger diagrams 



could be effectively reduced to smaller ones with fewer external lines. For the same result 
to occur in a general 0" theory, there must exist a general identity of the form shown in 
equation (100). We expect this identity to be true in all 0" theories, and we have explicitly 
computed the 0^ vertex 



Vabcd — 7 ^ 

{Aabc,d + ma + nib + m^ 



1 



njnblriclil - h + ma + Aa)„„^(l - h + mb + A6)_„j,(l - h + mc + Ac)-r. 

nid)m^—na—iii,—nc\^ab,c ~r ni^ -|- TJlb "Hlc fla Hb + Ti^jmc— ric— ria— ?i(, 



{nic -Ua-Ub- Ucy.nJnbK^c -Uc-Ua- nb)\ 

Vcdajma -ng-ria- na){Acda,b + mg - ng -Ug - rrih + Ub + rib) m^-n^-n.-ng 

{l-h + mb-nt, + Ab)-n,na\{l -h + mg-ng + Ag) _ria-nai^b -rib -rib -fig 
and checked our identity (100) numerically for the case of 0^ and 0^. 



We can use the functional equation to generate more general vertices. The idea is that 
if we construct the {s + t)-pt Mellin amplitude from [s + 1) and (t + l)-pt vertices and plug 
the result into the functional equation, then we will automatically generate an (s + t)-pt 
diagrammatic vertex. This means we begin by taking (schematically) 



M (6) = \^ VA-,...AsA{mi, ...,ms, m)SAim)VA^...AtA{ms+i 



,ms+t,rn) 



(167) 



where we are imagining that the external legs are themselves propagators coupling to further 
vertices, and 5 is the single propagator variable. Once we plug this formula into the functional 
equation, we need to deal appropriately with the fact that the external legs are 'off-shell'. 
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which means that the rrii can be shifted. The resulting identities can be used as a definition 
or a check of the n-pt vertex; we also obtain very non-trivial identities from the symmetry 
properties of these vertices. 



C.3 Flat Space Limit Identities 

In this appendix we prove two identities that were used in the study of the fiat space limit 
of the AdS factorization formula in section [6j Firstly, we want to prove that 



E 

m=0 



1 fd d Y {tLtaf-^-^e-^^/'^-Mtt 



m\ ydti dtji 



T{A-h + l + m) 



(3^L-\-f'^ 6{/3l - (3r) . (168) 



tT=tn=l 



This can be proven by integrating both sides against (3^ ^/3^ over (3^ and /3/j from to oo. 
After performing the integrals, the left hand side gives 



E 

m=0 

oo 

E 

m=0 



r(x)r(y) 



d d 



m!r(A-/i + 1 + m) \dtLdtR 
T{x)V{y) 



{tLtRf-'^-H^rty 



L'^R 



tL=tR = l 



m\T{A -h + l + m){x + h- A)_„(y + h- A)_, 



(169) 



=r{x + y + h - A - 1) 

and the right hand side trivially gives the same. 
Secondly, we want to prove that the limit 



lim > 

A-s>oo ^^-^ 



m=0 



d d 
dti dtR 



{PLPR)^{tLtR) 



h-A-1 -l3L/tL-l3R/tf, 



is given by 



{u + ^±i^)r(A -h + l + m)m\ 
2(3l 



tT=t„ = l 



5{Pl - PR)e-^'Pt 



2Plu + 1 



:i7o) 



(171) 



This was used in section 6.1.1 with the identification u = —6lr/A. 



The strategy to prove that ( 170 ) equals ( 171 ) is to compute their double Mellin transform 



by integrating them against (3^ /3|, over (3l and /3r from to oo. The Mellin transform 



of (171) is 



d/3e-^(3''+y+''-^—^ = / ds^—— ^—r- 

2l3u + l Jo s'^+y+hu + s/2 



(172) 
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where we changed integration variable to s = 1//3. The Melhn transform of (170) is 



hm > 

A-5-00 ^-^ 



m=0 

oo 



dtL dtR 



[u + 



A+2m' 



r(A -h + l +m)m\ 



(173) 



r(a: + t)r(^+f) 



tL=tR = l 

1 



liiii 7 



2m 
2 



Comparing with (172), we identify the continuous integration variable s with the limit of 



the discrete variable 4m/ A^. Then, the sum over m turns into the integral over s 



m=0 



A2 



— ^ / ds .. . 

'o 



(174) 



and, using the Stirling approximation, the summand reduces to the correct integrand. 



lim 



A2r(x + f)r(y + f) 



A^oo 4m!r(A - /i + 1 + m) (/i + x - f )_^ (/i + y - f ) 



-1/. 



sA2 



gx+y+h 



;i75) 
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